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Abstract 

Demanding a consistent quantum field theory description of spin ^ particles near a circu- 
lar Fermi surface in 2d leads to a unique fermionic theory with relevant quartic interactions 
which has an emergent Lorentz symmetry and automatically has an 5*^(4) = 5*0(5) in- 
ternal symmetry. The interacting theory has a low-energy interacting fixed point and is 
thus a non-Landau/Fermi liquid. Anti- ferromagnetic (AF) and superconducting (SC) or- 
der parameters are bilinears in the fields and form the 5-dimensional vector representation 
of 5*0(5). An AF phase occurs at low doping which terminates in a first order transition. 
We incorporate momentum dependent scattering of Cooper pairs near the Fermi surface 
to 1-loop and derive a new kind of SC gap equation beyond mean field with a d-wave gap 
solution. Taking into account the renormalization group (RG) scaling properties near the 
low energy fixed point, we calculate the complete phase diagram as a function of doping, 
which shows some universal geometric features. The d-wave SC dome terminates on the 
over-doped side at the fixed point of the RG, which is a quantum critical point. Optimal 
doping is estimated to occur just below 3/27r^. The critical temperature for SC at optimal 
doping is set mainly by the universal nodal Fermi velocity and lattice spacing, and is es- 
timated to average around \AQK for LSCO. The pseudogap energy scale is identified with 
the RG scale of the coupling. 
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I. INTRODUCTION 



In the renormalization group (RG) framework, Landau's theory of Fermi hquids 
is characterized by the irrelevance of the interactions of particles near the Fermi sur- 
face, in other words the low energy fixed point is simply a free theory of fermions. 
The underlying reasons for the wide success of Landau/Fermi liquid theory are well- 
understood , and consequently the known models of non-Landau/Fermi 
liquids are relatively rare. (Henceforth referred to simply as non- Fermi liquids.) An 
important exception is the Luttinger liquid and other related models consisting of 
quartic interactions of Dirac fields in d = 1 spatial dimension. Here the non-Fermi 
liquid behavior can be attributed to the fact that in Id, quartic interactions of Dirac 
fields are marginal operators in the RG sense. In higher dimensions quartic interac- 
tions of Dirac fermions are irrelevant and this is one of the reasons why candidate 
non-Fermi liquid models were not found in the previous works. Whereas more exotic 
non-Fermi liquid models have been proposed which typically involve gauge fields, 
the lack of non-Fermi liquid models in 2d appears paradoxical when one considers 
even the simplest models of itinerant electrons with quartic interactions, such as 
the Hubbard or t-J model, which are believed to be at strong coupling. Since such 
models have been proposed as good starting points for thinking about high su- 



perconductivity in the cuprates 



7|, |9|], it is certainly worthwhile to continue to try 
and construct relatively non-exotic models of continuum fermions with quartic in- 
teractions that have some resemblance to the Hubbard model and have non-Fermi 
liquid behavior in the normal state. 

Though the search for a novel kind of non-Fermi liquid in 2d provided one of the 
main initial motivations for the formulation of the model that will be presented and 
analyzed in this work, the model turns out to have many unexpected bonus features, 
almost all of which are intrinsic to 2d. We list the most prominent: 
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• The 4-fermion interaction is unique for spin | electrons and automatically has 
50(5) symmetry. In 2d the interactions are relevant and the model has a low en- 
ergy interacting fixed point with non-classical exponents which can be computed 
perturbatively. 

• The model generalizes to N fiavors, where N = 2 corresponds to spin ^ elec- 
trons, and has Sp{2N) symmetry. Since Sp{A) = 5*0(5) this provides a underlying 
framework based on a microscopic theory for exploring the ideas of Zhang based on 



50(5) IQillll- In particular one can derive the effective Ginzburg-Landau theory. 



• Because of the 50(5) symmetry the model naturally has both anti-ferromagnetic 
(AF) and superconducting (SC) order parameters that form the 5-dimensional vector 
representation of 50(5). For repulsive interactions the model has AF order and no 
SC order in mean field approximation. 

• When one incorporates momentum dependent scattering to 1-loop to go beyond 
mean field, an attractive d-wave channel opens up and the momentum dependence 
of the gap can be calculated. This d-wave SC phase terminates on the over-doped 
side at the RG fixed point, which is a quantum critical point. Due to mathematical 
properties of the d-wave gap equation, it also terminates on the under-doped side 
yielding a "dome". Due to the properties of the RG fiow, this attractive d-wave 
instability exists for arbitrarily strong repulsive interactions at short distances. 

• Although the model may be at arbitrarily strong coupling at short distances, the 
low energy fixed point is at a relatively small coupling 1/8, and this renders the 
model perturbatively calculable. We are thus able to calculate the main features of 
the complete phase diagram as a function of a doping variable, including the phase 
boundary of the d-wave superconducting dome and estimate the optimal doping frac- 
tion, which is near 3/27r^ ^ .15. This phase diagram depends on a single parameter 



< 7 < 1 which encodes the ratio of the strength of the couphng at short versus 
long distances. In the figure below, we summarize the results of our calculations for 
7 = 1, which corresponds to infinite coupling at short distances. The overall scale of 
temperature is set by the universal nodal Fermi velocity and the lattice spacing. 



Many of the above properties were highlighted on a list of the most iniportant 
features of high superconductivity compiled early on in the subject [g], IqI. Be- 
cause of the importance of the Cu02 planes, high is believed to be essentially 
a 2d phenomenon. This, along with the detailed properties of the solution of our 
model, in particular the phase diagram, led us to propose it as a model of high 



Tc superconductivity 121]. If our theory turns out to be the correct description, it 
reveals that the phenomenon of high superconductivity is remarkably universal, 
with a single energy scale, and its main features follow from the existence of the low 
energy fixed point in 2d. It is truly a beautiful phenomenon that has managed to 
realize some subtle theoretical loopholes in the usual requirements of unitarity, the 
spin-statistics theorem, and the Mermin- Wagner theorem, which are only possible 
in 2d. Our theory represents a significant departure from the models considered 
thus far in connection with high T^, which are reviewed in fisl . 14, Q, 16|, along 



with reviews of experimental results. On the other hand, we believe it represents a 
particular scaling limit of the Hubbard model at and just below half-filling, and is 
in this sense conservative in comparison with other more exotic ideas, and is thus 



in line with the ear 
Hubbard models 



sarly id( 

0,0,3. 



ideas concerning the role of AF order and the Heisenberg and 
However our model isn't simply a direct scaling limit of the 
Hubbard model with no attention paid to the Fermi surface, since the latter only has 
an S0{4) = SU(2) ® SU{2) symmetry, whereas our theory has the S0{5) symmetry. 
In our theory the "fermion sign" problem is solved by doing analytic, perturbative 
calculations in a fermionic theory from the beginning, and relatively simple 1-loop 




FIG. 1: Calculated phase diagram as a function of hole doping, which depends on a single 
parameter < 7 < 1 determined by the strength of the interaction at short distances. We 
set 7 = 1 corresponding to infinitely strong coupling. The vertical axis represents the low 
energy scale relative to the cut-off. What is shown are solutions 5'^ and 6'^ of the AF and 
d-wave gap equations (I109|120p in units of the cut-off Ac. The critical temperatures are 
proportional to these gaps with constants of order unity which we estimate, eqns. (I167|169]) . 
The overall temperature scale is determined by the universal nodal Fermi velocity and 
lattice spacing (jl71l) . The running coupling g in the gap equations is in terms of x in 
eq. (I147p . where x is the inverse dimensionless coupling. The hole-hoping h{x) is the one 
loop corrected expression (jl60p . The straight line Tpg is the renormalization group scale 
corresponding to the energy scale of the coupling, eq. (jl44p . and represents the boundary 
to the pseudogap region. The AF transition point at Haf = 3/47r^ is first order. The 
SC transition at /i* = 3/27r^ is second-order and corresponds to the fixed point of the 
renormalization group, i.e. a quantum criti^l point. The transition point hi ~ 0.13 is not 
universal, but relies on mathematical properties of the gap equation. 



calculations already reveal the main features. 

Irregardless of whether our model has been exactly realized in the laboratory, it 
can serve as a useful tool for exploring many of the paradigms in the area of strongly 
correlated electrons and also for developing new methods. For instance, we develop 
new gap equations that take into account higher order scattering of Cooper pairs near 
the Fermi surface. Our analysis shows clearly how in 2d one can obtain a momentum 
dependent gap with a d-wave structure from a rotationally invariant continuum field 
theory, i.e. without an explicit lattice that breaks the rotational symmetry. This is 
interesting especially since the precise origin of the d-wave symmetry of the SC gap 
has been unclear. We also show how to introduce doping in terms of the coupling 
and RG scale, and a small non-zero temperature as a relativistic mass coupling. 

For the remainder of this introduction we outline the organization of the paper 
and summarize our main results. In section II we motivate the model by showing 
how it can approximately describe particles and holes near a circular Fermi surface. 
The manner in which we expand around the Fermi surface is in the same spirit as 



in Q, y, y, y, e 



l5|] but differs in some important ways. For a single spin-less fermion 
one thereby obtains a free hamiltonian of particles and holes with a massless, i.e. 
relativistic dispersion relation. In section III we insist on a local quantum field de- 
scription of the effective theory near the Fermi surface with a consistent quantization. 
Since the particles are massless, the only known candidate field theories are either 
Dirac or "symplectic" fermions, which differ primarily by being first versus second 
order in space-time derivatives respectively. For the remainder of the paper we focus 
on symplectic fermions since unlike the Dirac fermions, the interactions are relevant. 



The model was first proposed in this context by one of us[17|], where the ground- 
work was done on the low energy non-Fermi liquid fixed point and in part the AF 
properties; at the time the SC properties were unknown. As explained in the present 
paper, the central idea of this previous work, that the AF order parameter is bilin- 



ear in symplectic fermion fields and that the low energy RG fixed point describes a 
quantum critical point, appears to be correct; however as we will see, the quantum 
critical point terminates the SC rather than AF phase. Quantum critical points in 



the context of high Tc were emphasized earlier by Vojta and Sachdev 19|]. The issue 



of the unitarity of our theory was mostly resolved in [18[ by noting that the hamil- 
tonian is pseudo-hermitian and this is sufficient for a unitary time evolution. In this 
paper, the expansion around the Fermi surface provides a new view on the pseudo- 
hermiticity and it is explained how it is related to the kinematics of particles versus 



holes. The critical exponents were computed to 2- loops in 18|, which corrected some 
errors in[l3]. In this paper we analyze many more properties, in particular the AF 
and d-wave SC ordering properties for the first time. 

Since the consistency of the quantization of a fermionic theory with a lagrangian 
that is second order in time derivatives is at the heart of the unitarity issue, in 
section IV we work out in detail the d = dimensional quantum mechanical case 
where all the subtle consistency issues are present. In this section we also construct 
the conserved charges for the Sp{2N) symmetry. In section V the field theory version 
in d spatial dimensions is defined and spin and charge are identified for the case of 
N = 2. In this section we also define the 5*0(5) order parameters for AF and SC 
order. 

In section VI we sketch an argument that the resistivity is linear in temperature 
in the limit of no interactions. In the next section we consider small thermal pertur- 
bations near T = 0. By comparing with the specific heat of a degenerate electron 
gas, we argue that a small non-zero temperature can be incorporated as a coupling 
in the lagrangian corresponding to relativistic mass m = aT and we estimate the 
constant a. 

In section VIII the mean field analysis is carried out with potential competition 
between AF and SC order. As we explain, these two phases actually do not compete 



in our model in this approximation. As a check of the formahsm, we reproduce some 
of the basic features of the BCS theory for an s-wave gap in the case of an attractive 
couphng in section IX. For repulsive interactions we find only AF order is possible in 
mean field approximation and this is studied in section X. There we first argue that 
this phase must be anti-ferromagnetic by comparing our model with the low energy 
non-linear sigma model description of the Heisenberg anti-ferromagnet. This gives 
another motivation for our model at half-filling away from the circular Fermi surface, 
and explains how the same model can interpolate between a SC phase and an AF 
one. We argue that the AF phase terminates in a first-order phase transition. The 
AF gap is then the solution of a transcendental equation that is analyzed in various 
limits. 

In section XI, orbital symmetries of momentum dependent gaps are studied in a 
model independent way and we explain how a d-wave gap can arise. This analysis is 
based on a gap equation which is derived in Appendix A. In section XII we compute 
the 1-loop contributions to the scattering of Cooper pairs and show that at low 
energies the d-wave channel is attractive if the number of components < 3. Since 
the theory is free for = 1, this means that only the physically relevant N = 2 
case has d-wave SC. This also means that the d-wave pairing cannot be studied with 
large A^ methods. 

Section XIII is devoted to describing our RG prescription which is specific to 2d. 
This is necessary for a proper understanding of the phase diagram. In section XIV 
we present global features of the phase diagram, which is characterized by some 
universal geometric relations, and bears a striking resemblance to the cuprates. The 
SC phase terminates at a second order phase transition precisely at the low energy 
RG fixed point, and is thus a quantum critical point. We also estimate optimal hole 
doping. In section XV we present detailed numerical solutions to the AF and SC 
d-wave gap equations at non-zero temperature. For reasonable values of the lattice 
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spacing and universal nodal Fermi velocity, we estimate Tc ~ 140/^ on average for 
SC in LSCO. In section XVI we describe the interpretation of the pseudogap within 
our model. 

Although we do not give a complete and rigorous derivation of our model from 
lattice fermion models, in order to motivate and point out relations, we have collected 
some known results about the latter in Appendix B. 

II. EXPANSION AROUND THE FERMI SURFACE 
A. Kinematics 

Let us first ignore spin and consider a single species of fermion described by the 
free hamiltonian in momentum space 



where /i is the chemical potential, we have defined (d^'k) = (i'^k/(27r)'^, and 



At finite density and zero temperature, all states with e < Sp are filled, where the 
Fermi energy ep depends on the density, and at zero temperature fi = ep- The Fermi 
surface Sp is the manifold of points kp satisfying e(kp) = ep. 

We wish to consider a band of energies near Sp a.^ shown in Figure [2] for d = 2. 
Let k be any wave- vector in such a band, and let r(k) denote a ray from the origin 
to infinity along the direction of k. We further assume that the Fermi surface is 
sufficiently smooth, such that r(k) intersects Sp only once. The latter implies that 
k can be uniquely expressed as 




(1) 



{4ck'} = (27r)'^5('^)(k-k') 



(2) 



k = k^(k) + p(k) 



(3) 
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where ki;'(k) is the vector from the origin to the intersection of r(k) with Sp- Whereas 
the two vectors k and kp'(k) are by construction parallel, the vector p(k) is either 
parallel or anti-parallel to k. Let us fix p to be a small vector parallel to k, i.e. point- 
ing radially outward, as shown in Figure O Since now ki?(k) is uniquely determined 
by p, we may write kp'(p). Furthermore, since the particles below the Fermi surface 
correspond to — p, the energies near the Fermi surface are approximately given by 

£(k) =£^±p-V^(k) (4) 

where ± corresponds to above or below the Fermi surface, and 

vp(k) = V£(k)|k, (5) 
is the Fermi velocity normal to Sp. 

Let us now assume that the Fermi surface is rotationally invariant, i.e. e(k) 
depends only on |k|. In 2d the Fermi surface Sp is thus a circle. This leads to the 
simplification that kp = \\<.p(k)\ and vp = |vi;'(k)| are independent of k. For any k 
in the band, 

k = {kp ± p)p (6) 

where p = |p|, and p = pp. Furthermore, since \p is normal to Sp, the energies are 
linear in |p|: 

£(k) = Ep ±vp IpI (7) 

For non-relativistic particles with e(k) = vp = kp/m^,. 

Since the map from k to p is one-to-one, we can define the following operators: 

Op = Ck^+p, bp = ci^„p (8) 

where it is implicit that kp depends on p. The above is a canonical transformation, 
since 

{alap,} = {fet,, V} = {27t'')6^'\p - p') (9) 
10 




FIG. 2: Expansion around the Fermi Surface 
After normal ordering, the hamiltonian for the particles in the band is defined to be 

H= (c/'^p) [{vfIpI - l2)al,ap + {vf\p\ + lippbp] (10) 

where = fi — ep is zero at zero temperature. For the remainder of this paper 
we mostly set /i = 0. Since we are only interested in a band of energies near the 
Fermi surface, we have introduced a cut-off A^. The vacuum |0) is defined to satisfy 
ap|0) = 6p|0) = 0. This corresponds to Ckj,_|_p|0) = c|^^_p|0) = 0, which correctly 
implies that all states with e < ep are filled. The Op and bp thus correspond to 
particles and holes respectively. 

There is an approximation made in obtaining the above hamiltonian having to do 
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with the density of states, and this is crucial to understanding how our expansion 
differs from previous works. In the rotationally invariant case, for particles above 
the Fermi surface: 



y"A = J dnj dp{kF + pY~^ [11] 



where dVL are angular integrals. Note that due to eq. (jH]), the angular integrals for 
k and p are identical. At least two approximations to the above are meaningful. 
The first favors low energies where one approximates k_p + p ~ kp- This is the 
approximation that is commonly made in the literature Q]. One the other hand, 
expanding out the {kp +pY~^, at high energies the leading term is p'^~^, and is the 
most sensitive to the short distance physics and spatial dimensionality. A possible 
short-coming of the first, low-energy approximation is that the high energy physics 
is discarded from the beginning. It cannot be recovered by the RG flow to low 
energies since the latter is irreversible. In the physical problem we are considering, 
the short-distance physics of the strong Coulomb repulsion is known to be important 
for understanding the AF phase, so it makes sense to adopt an approximation that 
favors high energies from the beginning and to then incorporate their effects by an 
RG flow to lower energies. We thus keep the most important term at short distances 
and set d'^'k = d'^p, i.e. kp + p ^ p. This is in line with the usual RG idea that it 
is important to fix the high-energy physics as accurately as possible, and then flow 
down to lower energies. Finally our choice is necessary for the 2d effective field theory 
description in the next section. However one should not conclude that every theory 
with a circular Fermi surface can be described by a relativistic field theory. One 
signature of a relativistic description is a density of states that is linear in energy: 
J d'^p = 2tivp J dee. 

The above expansion around the Fermi surface is thus not identical to the expan- 
sion m 

mm 

, where the integration over p is taken to be normal to Sp times 
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the angular integrations, and it is assumed that p <^ kp- This leads to the choice 
f (i'^k = f dVt f dpki^^, and the constant kp is absorbed into the definitions of the 

nnnnn 

operators. Thus in the approach followed in |2|, l3|, U, l5j, although the angular 
integrals obviously depend on d, the scaling analysis of the p dependence leads to 
marginal 4-fermion interactions for any d, and the resulting theory is effectively 1- 
dimensional, or a collection of such theories, one for each angular direction. Notably, 
it was not possible to obtain a non-Fermi liquid based on 4-fermion interactions in 
this approach {2]. 

In contrast, in the approach developed in this paper there is a strong dependence 
on d, as in other critical phenomena, and this will turn out to be very important. 
In particular, it leads to a non-Fermi liquid in 2d. There are other important justi- 
fications for this choice. In particular, near half-filling where the interacting lattice 
model can be mapped to the Heisenberg anti-ferromagnet, there is known to be a 
relativistic description of the low energy, long wavelength limit in terms of the 0(3) 
non-linear sigma model. It will be shown in section X that our choice of field the- 
ory near the circular Fermi surface can be extrapolated to half-filling in that an 
independent derivation of it can be provided exactly at half-filling. 

For general processes, physical momentum conservation of the k's is not equivalent 
to p conservation. However, consider a zero-momentum process proportional to 
kj). If the k's are all exactly on the Fermi surface, then this implies Yli ^fO^i) = 
0. For even numbers of particles, since all vectors on the Fermi surface have the same 
length, this is satisfied by pairs of particles with opposite kp- Allowing now small 
deviations pj from the Fermi surface, one has 

^(EkO=^( E P^-Ep^) (12) 

i particles holes 

Because of the particle/hole transformation for the 6's in eq. (IE]), this is equiva- 
lent to overall p conservation. Note that by construction, it is not possible for the 
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momentum of a particle and a hole to add up to zero, so in the above (5-function, 
holes are paired with other holes, and particles with other particles. Therefore, spa- 
tial translational invariance of our local field theory will ensure physical momentum 
conservation of the k's for this class of processes. 

The important allowed processes are shown in Figure [31 Examples of un-allowed 
process are shown in Figure HI The distinction between allowed and un-allowed pro- 
cesses can be made explicit by introducing an operator C that distinguishes particles 
and holes: 

CapC = ap, CbpC = -bp (13) 

where C is a unitary operator satisfying C = so that = 1. An eigenstate 
with pairs of particles and / or pairs of holes is then required to have C = 1. We will 
return to this in connection with the pseudo-hermiticity of symplectic fermions in 
the sequel. 




FIG. 3: Allowed processes 

B. Lattice fermions 

In the sequel, our field theory model will be related to lattice models of intinerant 
electrons such as the Hubbard model, although we do not claim a precise equivalence. 
The known 2d square lattice model results we will need to make the comparison are 
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FIG. 4: Un-allowed processes 

all contained in Appendix B. In this section we consider only the free, hopping term. 
In momentum space the 1-particle energy is fll9Up 

Ek = — 2t(cos k^a + cos kya) (14) 

where a is the lattice spacing. Equal energy contours in the first Brillouin zone are 
shown in Figure [5l The Fermi surface at half-filling is the square diamond with 
corners on the x, y axes. Note that one does not have to be very far below half-filling 
for the contours to be approximately circular. The free local field theory model in the 
next section can thus be viewed as an approximate effective theory for free particles 
on the lattice below half-filling. 

An important point is that our model is not simply a direct continuum limit of 
the lattice model since, without additional care, the latter does not take into account 
the Fermi surface at finite density. For instance, whereas the Hubbard model has at 



most an S0{4) = SU{2) SU{2) symmetry |36l]. our continuum model has the larger 
5*0(5) symmetry. Furthermore, as will be explained in section X, the success of our 
model can be attributed to the fact that an alternative justification of it can be given 
right at half-filling so that it can actually interpolate between half-filling and below. 
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FIG. 5: Fermi surface contours for free lattice fermions in 2d with lattice spacing a = 1. 
III. REQUIREMENTS ON THE FREE LOCAL FIELD THEORY 

The main requirements we impose for a local field theory description of the last 
section are: 

(i) The theory has a lagrangian description with a consistent quantization. 

(ii) In momentum space the hamiltonian reduces to equation fllOp for particles and 
holes of energy f^lpl- The latter is a relativistic dispersion relation for massless 
particles. 

In order to motivate our arguments, let us start from non-relativistic particles with 
= k^/2m^,. The second-quantized description consists of a single field \E'(x, t) with 
lagrangian 




(15) 



The field has the momentum space expansion 




(16) 
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Expanding around the Fermi surface one finds 

^(x, t) = e-''^' j (rf'^p) $p(x, t) e'^^^py^ (17) 

where 

^p(x, t) = Op e-^^^IPl*+^P-^ + bl e™^IPl*-^P-^ (18) 

We wish to find an effective theory for which satisfies the relativistic wave 
equation: 

(9,2-4v2)*p(x,t) = (19) 

Thus, due to the kinematics of the expansion around the Fermi surface we identify 
an emergent Lorentz symmetry. In 2d this Lorentz symmetry is 5*0(3). The case of 
d = 1 is special in that the Fermi surface consists of only two disconnected points 

= ±kp, and the decomposition (1171) naturally separates into left and right movers. 
For higher dimensions there is no such separation since all points on the Fermi 
surface are related by spatial rotations and continuously connected. There are many 
additional reasons why c? = 1 is the exceptional case, and these will be pointed 
out where appropriate in the sequel since some of the literature attempts to draw 
analogies between Id and 2d. 

There are only two known candidate field theories which differ in whether the 
lagrangian is first or second order in derivatives. First consider the case of first- 
order. One then needs to factor the operator — into two first-order multiples. 
The only way to accomplish this is to promote \E' to a multi-component field and 
introduce a matrix representation 7^, yU = 0, 1, .., d, of the Clifford algebra: 

{7^7^ = 2//^^ (20) 
where r]'^'^ = diag(l, —1, ., —1). One then has 

J2 1^1'' 9 A = E = - V' (21) 
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(We have adopted the relativistic notation = {x^, x^, .., x'^) = {t, x) and 9^ = 
and henceforth, repeated indices are imphcitly summed over.) The lagrangian is then 
the standard first-order Dirac lagrangian: 

d'^:>^i^Ydf,iJ (22) 

where ip = ip'^'-f'^. The smallest representation of the Clifford algebra is two dimen- 
sional: 7° = cr^,7^ = iCx, and 7^ = iay where a are the standard Pauli matrices. In 
this simplest case, although there is a doubling of components, they are constrained 
by the Dirac equation of motion and the spectrum still consists of one species of 
particles and holes with hamiltonian eq. ( ITOl) . 

Interactions are relevant to the low energy physics if the operator characterizing 
them has scaling dimension less than d + 1. The classical scaling dimension of the 
Dirac field ip is d/2 in d spatial dimensions, thus a quartic interaction has dimension 
2 ■ d and is thus only relevant for d < 1. It is not even perturbatively renormalizable 
for d > 1. Thus, the Dirac theory should lead to ordinary Landau- Fermi liquid 
behavior in 2d. It is noteworthy that once again (i = 1 is special, and this helps to 
explain how for instance the Id Hubbard model can be mapped onto interacting Dirac 
fermions, and the low energy fixed point found using special bosonization techniques 



and spin-charge separation 35|. Simply based on the fact that the interactions in 
the Hubbard model are strong in 2d, one can rule out a description in terms of 
Dirac fields with quartic interactions since the latter are irrelevant. Furthermore, it 
is already understood that one normally needs additional special properties in order 
to obtain the first-order Dirac theory. For example, it is known to arise when one 
expands around special nodes (Dirac points) on the Fermi surface for a hexagonal 



lattice [2l|] , as in graphene, and the multiple components of the Dirac field correspond 
to different sub-lattices. There is no reason to expect this here for a square lattice, 
as in the cuprates. 
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The other candidate field theory is second-order in derivatives, with kinetic term 

S = [ dtd''^d,x-d^X^ (23) 



For fermionic (Grassman) fields, this is a very unconventional theory, since it po- 
tentially has problems with the spin-statistics theorem and unitarity; in high energy 
elementary particle theory it usually corresponds to ghost fields. These issues will 
be discussed in detail and resolved completely in the next two sections. Here, let us 
give the main arguments for why this should be the right starting point: 

(i) As shown in the next two sections, the free theory in momentum space corresponds 
precisely to the hamiltonian ( fTOl) for particles and holes near a circular Fermi surface. 
This is of course a perfectly hermitian theory with no negative norm states. 

(ii) The fundamental field x has scaling dimension {d — l)/2 and thus quartic inter- 
actions have dimension 2{d — 1) which is actually relevant for d < 3. Thus it can 
have non-Fermi liquid behavior. 

(iii) Although we are led to consider this model for the nearly circular Fermi surface 
below half-filling, a simple argument leads to the same model at half-filling. It is 
well-known that a low-energy description of excitations above the staggered AF state 
is described by the 0(3) non-linear sigma model for a field (p constrained to have 
constant length with the action 

S = [ dtd'^xdj-d''$ (24) 



(See Appendix B.) In our model the anti- ferromagnetic order parameter (j) is bilinear 
in the fields = X '^X^/v^- The non-linear constraint on the fields follow from 
imposing a similar constraint on the x fields: X~X~^ = constant. This was pointed 



out in[l7j. Inserting this into the above action one finds that one obtains the second 
order action fl23|) for the x fields up to irrelevant operators (eq. (11041) below). Thus 
the symplectic fermion model with interactions can in principle describe AF order, 
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and in the sequel we will show that this is indeed the case. This is explained in 
greater detail in section X. 

Since vp only serves to convert dimensions of time and space, it plays the role 
of the speed of light; we can set it to unity since it can always be restored by 
dimensional analysis. The important point here is that vp is fixed and universal in 
our model, in particular it does not depend on the coupling. In the sequel the doping 
will be related to the coupling, so vp does not depend on doping either. For high 
Tc materials, we believe the "speed of light" has actually already been measured 20| . 
The Fermi surface at half-filling is closest to the nearly circular surface below it in 
the nodal k = (0, 0) to (vr, tt) direction, so our vp should correspond to this nodal 
Fermi velocity. Remarkably, the latter was measured to be universal at low energies. 



i.e. independent of doping in 



201 ]. Whereas this universality oi vp has not been 



explained theoretically up to now, it is a necessary aspect of our theory. Taking the 
slope dE/dk of the curves ir^^ we estimate vp ~ lAevA = 210km/ s for LSCO. 
As we will show in section XV, this gives very reasonable estimates of Tc. Another 
signature that the system may be in a relativistic regime is a density of states that 
is linear in energy, as explained in the last section. 



IV. SYMPLECTIC FERMION QUANTUM MECHANICS 

As stated above, symplectic fermions are primarily characterized by a lagrangian 
that is second-order in space and time derivatives. Since this is unfamiliar to most 
readers and there are some delicate issues in the quantization of such theories, let us 
first start with the simplest case of d = quantum mechanics. 
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A. Canonical quantization 



In order to draw comparisons, let us first consider a first order lagrangian as in 
the Dirac theory: 

TV 

L = ^ {icldtCc, - ujclca) (25) 

a=l 

It is well-understood that this lagrangian has two consistent quantizations, i.e. 
one can impose either canonical commutation relations [cq, cj^] = 6a,i3 or canoni- 
cal anti-commutation relations {cq, cj^} = Sa^p- In both cases the hamiltonian is 
H = J2a^^a^a- K's clear that both options are possible since eq. ffT^ is a proper 
second quantized description of either bosons or fermions. For future reference we 
note that the model has a manifest SU (N) symmetry. 

The second order bosonic version of the above is just the ordinary harmonic oscil- 
lator with L = {{dtqY " t^^)/2. Since the first-order lagrangian can be consistently 
quantized as a fermion or boson, one expects that the second-order case should also 
be quantizable as a fermion, and as we now describe, this is indeed the case. In order 
to have a fermionic version, we need at least 2 degrees of freedom since fermionic 
variables square to zero. Let us therefore consider the lagrangian 

L = X~X'^ - ^'^X'X'^ (26) 
where x s-'^s Grassman variables: 

{X\X'} = 0, (27) 

which implies (x~)^ = (x^)^ = 0) ^-i^d we have defined x = 9tX- The canonical 
momenta are p~ = dL/dx~ = X^ and p"*" = dL/dx^ = — x"? which leads to the 
canonical anti-commutation relations 

{x',X^} = -{x^,r} = ^ (28) 
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The canonical hamiltonian is simply 

H = x'-t + ^\-X'^ (29) 

The equation of motion is [d^ + = 0. Because this is second-order, the mode 
expansion involves both positive and negative frequencies: 



= -L (-6^ e-'"* + a e'"*) (30) 



The canonical anti-commutation relations (l28l) then require 

{a,at} = {5^5t} = 1 (31) 

with all other anti-commutators equal to zero. The hamiltonian is 

E = uj{a^a + b'^b - 1) (32) 

B. Pseudo-hermiticity 

The only subtle aspect of the above quantization is the extra minus sign in the 
expansion of sq. (!30l) . which was necessary in order to have the canonical 

relations fl3T]) . This minus sign implies that x^ is not the hermitian conjugate of 
X~ ■ One can understand this feature more clearly, and also keep track of it, with the 
operator C that distinguishes particles and holes in eq. ( fT3l) : 



X 



C{x-yC. (33) 



In terms of the original x variables, the hamiltonian is pseudo- hermitian, = CHC. 
However after using the equations of motion and expressing it in terms of a, 6's, since 
it is quadratic in 6's, the hamiltonian ( !32l) is actually hermitian. This issue will be 
revisited when interactions are introduced in the next section. 
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C. Symmetries 



We now study the symmetries of the A^-copy theory. Introduce variables x^, 
i = —, +, a = 1,2, .., N and define the lagrangian 



L = lj2'^^i^^^--'''^-Xi) (34) 
where e^j is the 2x2 anti-symmetric matrix e |_ = — e+_ = 1. The hamiltonian is 

H=lj2'^draXi+^\aXi) (35) 

Arrange Xa i^^o ^ 2A^-component vector and consider the transformation x ~^ Mx, 
where M is a 2A^ dimensional matrix. Then the lagrangian is invariant if M^ej^M = 
e^v where = e® Iat, and M* is the transpose. This implies that M is an element of 
the group Sp{2N) of dimension A^(2A^ + 1). Interestingly, of the classical Lie groups, 
Sp{2N) is the only one that doesn't play any known role in elementary particle 
physics 251]. Note that the bosonic version of the theory with 2A^ real components 
has the symmetry 0{2N) with dimension N{2N — 1) + 1, thus the fermionic version 
always has a larger symmetry. The N-component symplectic fermion also has a 
larger symmetry than the N-component first order fermionic action which has a 
U{N) symmetry or an 0{2N) symmetry if the complex fermions are rewritten in 
terms of 2A^ real fields. 

The conserved charges that generate the N[2N + 1) dimensional Lie algebra of 
Sp{2N) are easily constructed. Define 



Qal3 = -^{XaX]} +X]}Xa) (36) 
Qal3 ~ ^ iXaXp + XpXtt) 
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One can easily verify that the hamiltonian H commutes with all the charges Q using 

{Xa^Xp = -{xt,Xp} = ^^afS 

{xl.,X^} = {XL,X^^} = (37) 
In "momentum" space, the charges are 



Qal3 


= aia^ - h^pba 




Qa/3 


= a)j)p + a^phc, 


(38) 


Qal3 


= b^ap + h^pa^ 





and they satisfy the hermiticity properties 

{Q%)^ = Ql, (Qap)^ = Qlp (39) 

V. FIELD THEORY VERSION WITH SPIN AND INTERACTIONS 
A. Lagrangian and Hamiltonian 

The field theory version in d spatial dimensions follows straightforwardly from the 
above d = case with the addition of spatial or momentum integrals. Introducing 
fields Xai^yt): the action is 

S=]^jdtd''^Y. {d'Xad.xi - m^xLxy (40) 

and the equations of motion are 

{d^d^ + m^)x = 0. (41) 
The momentum space expansion is 

X+(x, t) = I {-bl e-^^- + ap e^^-) (42) 
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where Up = a/p^ + m'^ and p ■ x = ujpt — p ■ x. (We do not display the a indices 
since they just correspond to identical copies.) The canonical anti-commutation 
relations are 



{x-(x,t),x+(x',t)} = -{x+(x,t),r(x',t)} = z5W(x-x') (43) 
which in momentum space leads to 

{ap,al} = {bp,bl'} = i^^rS^'\p-p) (44) 

(All other anti-commutators are zero.) 
The hamiltonian is 

N 



H = (d'^p) ujp (aj,,„ap,„ + 5{,,„6p,c,) 



(45) 



In the limit m ^ 0, we obtain the effective hamiltonian near the Fermi surface in eq. 
(ITOl) . as desired. The mass m in this section is an infra-red regulator and is unrelated 
to the non-relativistic mass m^, above. In section VII we will show that it can be 
viewed as proportional to the temperature. 

The field theory has the same Sp{2N) symmetry as the quantum mechanical 
version. The expressions for the conserved charges are identical to the eqs. fl36|38p 
with additional integrals over x or p. 

There is a unique 4-fermion interaction that preserves the Sp{2N) symmetry: 

^int = -n'g I dtd^^ {x-eNX^Y = -^T^^g I dtd"^ (^^^^^) ^^^^ 
(We have included an overall tt^ so that the RG equations below for q have no vr's; 

n 

our convention is the same as in [18|.) For N = 2, even without imposing the 
5*^(4) symmetry, there is a unique interaction due to fermionic statistics since there 
are only 4 independent fields, which implies that higher order terms beyond quartic 
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interactions are zero by Fermi statistics. This interaction is automatically 5*0(5) 
invariant. Positive g corresponds to repulsive interactions. Since the field x has 
classical scaling dimension [d — l)/2, the interaction is a dimension 2{d — 1) operator 
which is relevant for (i < 3 as previously mentioned. 



B. Pseudo-hermiticity 

The symplectic fermion action ( HOl) has a Lorentz invariance if x is understood to 
be a Lorentz scalar. Since x is fermionic, to a particle physicist this model would 
appear to violate the spin-statistics theorem. There are two separate aspects of this 
issue. First of all, in the condensed-matter context, rotational spin is an internal 
fiavor symmetry (spin | fermions corresponds to = 2) which is not viewed as 
embedded in the Lorentz group. This implies that spin | particles are not forced 
to be described by the first-order Dirac theory. There is no violation of the spin- 
statistics connection in our theory since we are quantizing spin ^ particles with 
fermionic fields. The potential problem rather has to do with unitarity, as explained 
in the quantum mechanical case studied in the last section, and is manifested in the 
pseudo-hermiticity property ( |33l) . This explains how the proof of the spin-statistics 
connection is circumvented: the proof assumes that the hamiltonian is built out 
of fields and their hermitian coniugates and thus doesn't allow for different fields 

n 

being related by pseudo-hermitian conjugation[5|. Furthermore, in the end the free 
hamiltonian in momentum space is a perfectly hermitian theory with no negative 
norm states. 

Whereas the free theory is hermitian in momentum space, for the interacting 
theory, it follows from ( l33l) that the hamiltonian is pseudo-hermitian: 

= CHC (47) 
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where the unitary operator C = and = 1. This sort of generahza- 
tion of hermiticity was understood to give a consistent quantum mechanics long 



ago by Pauh 22 L and more recently in connection with VT symmetric quantum 
23l . l24l |. Let us summarize the main properties enjoyed by pseudo- 



mechanics 
hermitian hamiltonians: 

(i) Define a C- hermitian conjugation as follows: 

A^^ = CA^C (48) 

Then the usual rules are satisfied: 

{AB)^^ = B^^A^% {aA + hB)^^ = a*A^^ + b*B^^ (49) 

(ii) Define a C-conjugate inner product: 

(50) 

Then time-evolution is unitary: 

{i^'{t)\m)c = (V^'|e^^^*Ce-^^*|^) = (^'(0)1^(0)), (51) 

(iii) The eigenvalues of H are real: 

{E - E*){ijE\i^E)c = {^Pe\{CH - H^C)\tPE) = (52) 

(iv) Diagonal matrix elements of C-pseudo-hermitian operators A = J^" are real. 
This will be important in the sequel since it guarantees the reality of vacuum expec- 
tation values of pseudo-hermitian order parameters. For our model 

H'^^ = H, (x")t^ = x+, (x^X+)^= = (53) 
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In the present work there is a new aspect of the pseudo-hermiticity that relates to 
the kinematics of the expansion around the Fermi surface. As explained in section II, 
for the 4-particle processes near the Fermi surface, conservation of p relative to the 
Fermi surface is equivalent to conservation of the physical k momentum if particles 
are paired with particles and holes with holes. Since C = ±1 for particles versus 
holes, on physical grounds we should restrict to eigenstates with even numbers of 
holes and even numbers of particles with C = 1. Let {iPe) denote an eigenstate of H 
which is also an eigenstate of C. Then H^iPe) = C^H\iIje) = H\iIje)- Thus for the 
eigenstates of interest, H = W . 



C. Charge and Spin 

In order to describe spin | electrons, as usual we treat the spin as a flavor and 
thus consider the N = 2 theory. The symmetry of the free theory is 5*^(4) = SO (5). 
A subgroup of this large (10-dimensional) symmetry can be identified with rotational 
spin and charge. 



It was pointed out in [18l] that there are potentially two ways to identify electronic 
spin, and the focus in that work was the SU{2) subalgebra that exists for all and 
acts on the ± indices of xt- turns out that the other identification is more natural 
in the present context since, as explained in section II, the two components corre- 
sponding to the ± indices are already necessary for the expansion of a single spinless 
fermion near the Fermi surface. There is actually an analog of the identification of 
spin and charge for arbitrary A^. Let M = e~"^ be an element of the group Sp{2N) 
in the defining 2N dimensional representation. Using the relation 

M^€nM = eN, (54) 
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the elements of the Lie algebra satisfy 

m*eN = — CArm, (55) 

and a basis is the following: m G {1 ® a, a^,. ® a^, ® a^;, o"z ® a^}, where cxj are the 
Pauli matrices, a is an x dimensional anti-symmetric matrix, and are N x N 
symmetric matrices 25|] . Clearly Sp{2N) has an 5'f/(2)®^ sub-algebra generated by 



ai (g) where J(°) = diag(0, 0, .., 1, 0, .., 0) with the 1 in the a-th entry. For = 2 
this corresponds to an 50(4) = SU{2) ® SU{2) symmetry. However since this sub- 
algebra does not mix the flavors, it is not the right sub-group for identification with 
spin as a flavor. The correct identification is the SU{N) sub-algebra generated by 
1 ® a and az^az where now is also traceless. There is also a f/(l) which commutes 
with the SU{N) corresponding to (Jz ® Ia?. The SU{N) is generated by the charges 
which form a closed algebra: 

The U{1) we identify with electric charge is generated by Qe = J2aQaa- 

Let us now return to the physically interesting case of = 2 and label the 
two components as a =], |, corresponding to up and down spins. The SU{2) spin 
symmetry is generated by 



= I Wt - Qli) = -f y ixjxt + xUj - xjxt - xUj) 

= |(ci'^p)(l(a^^apt-a;,^api)-i(6^^6pT-6^^6pi)) (57) 
= 73^n = -73 / ix^^t + XW = ^ / (d'pH.a,, - bib,,) 

satisfying the SU (2) Lie algebra: 

[Qz.Q±] = ±Q±, [Q+,Q-] = Qz (58) 
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(As before x = 9tX-) The U{1) charge is generated by 

Qe = g?t + Qll = -'J (xfXT+ + Xfxf + Xlxt + XUI) 

= jid'p) ((a^^apt + a;,^api)-(6;,^6pt + 6;,^6pi)) (59) 
The conserved electric current corresponding to the above charge is 

The fields have electric charge Qe = ±1. Commutations of the fields with 
the SU{2) generators shows that (xfjXj^) form a doublet whereas ix^^X'l) is the 
conjugate: 

[Qz. xf ] = T|xf , [Qz, xf ] = ±|xf (61) 

Note that the above identification is consistent with a being particles and h holes: 
Qz\ci'\) = Ikt) whereas Qz\b]) = — Also, the a-particles have Qe = I whereas 
the b have Qe = —1. 

The additional 6 conserved charges with electric charge ±2 that complete the 
50(5) Lie algebra are 

g± = -2t J d^'^xfxf, Qfi = -2^ / d'^xfxf (62) 
Qf, = -^j d'^ {xfxt + Xfxt) (63) 
These symmetries fiip the charge and spin of the fields, for instance: 

[Qtv ] = 2xf , [Qtv ] = xt. [Qtv xj] = xt (64) 

There is no real separation of spin and charge degrees of freedom in our model at 
this stage, unlike Dirac fermions in Id, where this separation relies on bosonization. 
However as we explain in section VIII, the Goldstone bosons are electrically neutral 
but carry spin quantum numbers. 

30 



D. SO {5) order parameters 



The 4 fields xf^ xf the 4-dimensional spinor representation of SO (5). (See 



e.g. 



25|). The important operators are bihnears, which decompose as 4 ® 4 



1 © 5 © 10, where 1 is the singlet x~X~^j the 5 is the vector representation, and 10 
the adjoint. The 10 corresponds to the currents constructed above. The 5 will serve 
as the order parameters of our model in the sequel. A triplet of fields (j) = {(pxi 'Py, 4'z) 
transforming under the 3-dimensional vector representation of the spin SU{2) is the 
following: 

0+ = XfXi+, <P' = Xlxt, = 71 (xf - Xrxf) (65) 

where (j)"^ = {(px ± i<Py) / V^- Note that these fields are electrically neutral. Let us 
also define two SU (2) singlets that carry electric charge Qe = ±2: 

(Pt = xtxt^ C = XfXf (66) 

Using the commutation relations of the 50(5) charges with the x-fields, one finds 
that the 5 order parameters 

$ = (0,,0y,0„0+,0;) (67) 

transform under the 5-dimensional vector representation of SO (5). For instance 
[Q-\P 4'e] ~ and [Q^p 4>~^~\ = 20^. The ordering of fermionic operators was cho- 
sen such that {^Y" = $, more precisely (p"^" = cj) and (0^)^'' = 0", which guarantees 
that vacuum expectation values of are real and those of 0^ are complex conjugates. 
An 50(5) invariant is the following 

$ . $ = 0+0- + 0-0+ + 02 _ 0+0; _ 0^0+ (68) 
Finally the interaction lagrangian density eq. 0461) can be expressed as 

Ant = ^ $ ■ $ = -87r2(^ X^xUat (69) 
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E. Low energy fixed point 



The Feynman rules for the theory are the same as for a Lorentz invariant scalar 
with 0^ interaction, which in practice is quite simple [s], except for some all 
important minus signs. In the sequel we mainly work in euclidean space t —it since 
this simplifies solving the gap equations and is also appropriate for finite temperature. 
Where appropriate we will return to real time (Minkowski space) for certain physical 
quantities. In particular the propagators in euclidean space are 

{xMxtm = -ixl{x)x-,m = -S., J (0je--- (70) 

where D = d + 1, and they respect causality. In the usual euclidean conventions the 
4-vertex for the interaction in ( l46l) is —An'^g. (See Figure O) 

Our RG prescription will be specialized precisely to 2d and is described in detail in 
section XIII after some 1-loop diagrams are explicitly calculated. In order to obtain 
a clear physical picture, it will turn out to be important to carry out the RG directly 
in 2d. However the RG can also be studied perturbatively in an epsilon expansion 



around d = 3. For future reference we summarize the main results obtained in |17l.ll8[| 
from the epsilon expansion. We will not use these results very much in the sequel; 
they are included as a guide to the sign and strength of the anomalous corrections 
to scaling. Since the perturbative expansion differs from that of bosonic scalars 
only by fermionic minus signs, the standard methods for the bosonic 0(M) vector 
models apply. Specifically, these are models of an M-vector of scalar fields (p with 
interaction (^ ■ 0)^. In D = d + 1 = 3, the low energy theory is the Wilson-Fisher 



fixed point [271] describing classical 3d magnets. The anomalous dimensions of the 
most important operators were computed to 2-loops in ISj. As it turns out, most 
of the results can be obtained from known results for the 0{M) models analytically 
continued to M = —2N. Besides the x fi^ld itself and mass term x~X'^ (which 
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will be the thermal perturbation), an additional important class of operators are 
the bilinears that correspond to order parameters, in particular the cf) of the 5*0 (5) 
order parameters when = 2, of the form x~(^X^ where a is a Pauli matrix. These 
have no analog that is known to be physically meaningful for the 0{M) models. 
It should be emphasized that the quantum critical exponents of our N = 2 model 
are completely unrelated to the usual Wilson-Fisher exponents of the bosonic 0(3) 
model, since here the 0(3) vector is a composite bilinear operator and it's scaling 
dimension must be calculated in the fundamental x-theory. 

Let |X] denote the scaling dimension of X in inverse length, i.e. energy units. 
One can also define a correlation length exponent based on the mass m: 

e ~ rn-" (71) 

where u = 1/ [m] : 

u~' = (d+l-lx-xn)/2 (72) 



Specializing the results in 18| to = 2 and 2d one obtains 



1 ^ ^ 1 

M^^, lx~ffxn^^. ^^Y^ (73) 

Note that whereas x decreases in dimension from the classical value 1, x~^X~^ 
increases. 

Some remarks concerning the d = 1 case are again appropriate. It is well-known 
that the low energy fixed points of the bosonic 0(M) vector models do not extend 
down to d = 1, i.e. D = 2. For general M this can be viewed as a manifestation 



of the Mermin- Wagner result (29|, which states that spontaneous symmetry breaking 
in a D-dimensional euclidean field theory is not possible in D = 2. On the other 
hand, the models do have a conformally invariant fixed point for —2 < M < 2 in 



2D 



281 . Since our A^-component symplectic fermion model is formally equivalent 



in perturbation theory to 0(— 2A^), this suggests that our models do not extend to 
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d = 1 except possibly for — 1 < < 1. This helps to explain why for instance it has 
never been considered before in connection with lattice fermion models in Id. 
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VI. RESISTIVITY IN THE NORMAL STATE 



In this section we give a rough calculation of the temperature dependence of the 
conductivity when the interactions are negligible. The rigorous study of this question 
requires a full finite temperature treatment of the Kubo formula, which is known to 
be quite subtle, and therefore beyond the scope of the present article. We hope 
to return to this in a future work, however in this section we present the following 
non-rigorous scaling argument. 

Start with a version of the zero temperature Kubo formula 

a,,{q,u;) = - H dte''^' [ d^^e-'^'^ ([J,(x, t), J,(0)]) (74) 
^ Jo J 

where i,j = x,y are the spatial components of the currents given in fl60l) . 

It is well known that the frequency a; = 0, which corresponds to the DC con- 
ductivity, is a delicate limit. Let us first set u; = in the integrand. Using the 
propagators in eq. flTOj) one finds 

where we have set m = 0. The above integral is proportional to A^. Since the 
conductivity cr is dimensionless in 2d, one must have that it is proportional to Ac/T 
where T is the temperature. For reasons that are unclear, this apparently amounts 
to setting the overall in eq. ([71]) equal to T. Taking into account the two spin 
components, and the 4 terms of the above form, one obtains 



4 



Therefore the resistivity (l/o") is linear in the temperature T. 
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VII. THERMAL PERTURBATIONS AND ANOMALOUS SPECIFIC 
HEAT 



The parameters of our model thus far are the Fermi velocity vp, which we have 
set equal to 1, the cut-off A^, the coupling g, and the infra-red regulator mass m. 
The Fermi energy ep can be viewed as implicit in the cut-off if the latter is taken to 
be the frequency ujd = vphp. The density of the free electron gas is only a function 
of kp: 

Thus the density can be varied by varying the cut-off, and since g is proportional to 
the cut-off in 2d, equivalently by varying g. (See the discussion of the RG in section 
XIII.) 

In this section we suggest how to introduce a small non-zero temperature as the 
mass m. This is similar to how temperature appears in the Landau theory for con- 



tinuous phase transitions in 0{M) magnets |30|, where there also temperature corre- 
sponds to a coupling in an effective action. However there are some important differ- 
ences, since in the latter the coupling is T — T^, whereas here it will be proportional 
to T^. Since we have built our model by expanding around the zero temperature 
Fermi surface, it is not obviously consistent to incorporate a finite temperature by 
starting over and formulating a finite-temperature version of the x fields with Mat- 
subara frequency summations, ignoring that they are effective fields. Rather, since 
a small non-zero temperature amounts to a small distortion of the Fermi surface, it 
could correspond to an additional coupling in the lagrangian. If this is correct, then 
it should be possible to obtain known results in the limit where the interaction is 
turned off. Suggestive of this possibility is the fact that in natural units the leading 
contribution to the specific heat of a degenerate electron gas can be expressed entirely 
in terms of kp and T, i.e. the dependence on the electron mass is only through the 
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Fermi velocity vp = kp/m^. We emphasize that this way of introducing temperature 
is expected to be vahd for temperatures near zero, and cannot replace a full-fledged 
finite temperature formalism at arbitrary temperatures. It will however be useful for 
exploring the temperature dependence of low temperature gaps in the sequel. 

In euclidean space, introducing a finite temperature T is known to correspond to 
compactifying the euclidean time to a circle of circumference 1/T. Let us consider 
adding a thermal perturbation to the euclidean action of the form 



SS = j dtd'^x grOrix) (7^ 



where qt is a coupling and Ot is the "thermal operator" . Since the euclidean func- 
tional integral is over e~^, the correction to the free energy F = —TlogZ to lowest 
order in g-p is 6F = Tgx J dt(fi^{OT)- Since the euclidean space-time volume is 
Y^^) = V/T where V is the spatial volume, one finds 

^ = 9T{Or) (79) 
This leads us to identify the thermal perturbation with a mass term: 

gr = {aT)\ Ot = X^X^ = J2 ^« (^0) 

a 

The above relation is consistent with dimensional analysis for a a dimensionless 
parameter. The specific heat Cy at constant volume is Cy = — T|^. Since the 
propagator goes as — as m — (See eq. (jTOl) above) one has 

d^p 1 



iOr) - EixZx:) - -2 1 p^^. (81) 



where the factor of 2 comes from spin up and down. If the cut-off Ac is equated with 
the frequency kp, then this leads to specific heat that is linear in T: 

7y _ a^T kp 



Cv _ a^T kj-^ 

~V ~ {d - l)2'^-%('^+i)/2r(^^) ~p 



(82) 
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(We have re- introduced the Fermi velocity for the sake of comparison.) The above 
dependence on T, kp and vp is the correct one, i.e. it is the same as for a non- 



relativistic degenerate electron gas near T = 0. Repeating the standard calcu 



ation 



of the specific heat of the electron gas to order T/ep (see for instance |30|]) for 
arbitrary d, expressing the result in terms of kp instead of the density, and requiring 
the result to match eq. (182!) fixes the constant a: 

^^^.V-i)r(^) ^^^^ 

For d = 2, a = TT V6 ^ 1.7. The case (i = 1 is noteworthy since a = 0, and will 
be commented on below. 

The kp dependence of the specific heat in eq. f l5^ is a direct consequence of the 
scaling dimension d — 1 of Op, which is twice the scaling dimension of x- When one 
includes the quartic interaction, at the low energy fixed point this scaling dimension 
has anomalous corrections and this should lead to anomalous T dependence of the 
specific heat. The scaling dimension of Op was computed perburbatively in an 
epsilon-expansion results summarized in the last section for arbitrary 

A^. For N = 2, Op has scaling dimension approximately equal to 5/8 in the epsilon 
expansion. 

We can provide a naive estimate of the anomalous T dependence of the specific 
heat. Let {Op} denote the scaling dimension of Op. Then if we assume the only 
effect of the anomalous corrections is to replace k'j^T^ by k'^^'^\ then since Cy/V has 
scaling dimension d, this requires Cv oc T'^-I^rl. It is not clear that our previous 
assumption is correct however, since in general Cy could contain terms m^k^p with 
x + y = d. Nevertheless, using our estimate of 5/8 for the dimension of Op, this gives 
Cv oc T^^/^. Thus Cv/T oc T^/^. Though the exponent 3/8 should perhaps not be 
taken as very accurate, the general point is that since {Op} is shifted downward to 
5/8 from the classical value 1, this shows that Cy/T should vanish as T ^ 0. This 
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shift downward is entirely due to a fermionic minus sign 171]. We hope to study this 
more carefully in future work. 

Finally it is important to note that the manner in which we have introduced 
temperature allows, at least computationally, for phase transitions that break the 
continuous SU{2) and U{1) symmetries in d = 2. The Mermin- Wagner result 29| is 
the statement that spontaneous symmetry breaking is not possible in a 1 + 1 dimen- 
sional (space plus time) quantum mechanical system at zero temperature because 
of infra-red divergences that plague the existence of Goldstone bosons. An example 
of such a divergence is in eq. ( IHTi) for = 2. In the Matsubara approach to finite 
temperature, time is compactified into a circle of circumference 1/T and the discrete 
Matsubara frequencies are summed over. Thus the arguments of the theorem in 
principle apply to a finite temperature system in 2 + 1 dimensions at each Matsub- 
ara frequency. On the other hand, in our model temperature appears as a coupling 
in the theory, as in classical statistical mechanics in 3d. Furthermore, as discussed 
in section V, our model breaks down in d = 1, since the fixed point is lost, and a 
manifestation of this is the vanishing of a in d = 1. 

In the sequel, wherever m is non-zero, it should be thought of as representing a 
small non-zero temperature. 



VIII. MEAN FIELD ANALYSIS 

Because of the 5*0(5) symmetry, the interaction term in the lagrangian can be 
expressed in terms of either the magnetic or electric order parameters: 

/^int. = -8n^g X^xtxjxt = 8^'9 ■ 0/3 = Sir^g <Pt<P~e (§4) 

This implies that magnetic and SC order may in principle compete. In this section 
we study this in mean field approximation. 
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Introduce auxiliary fields s, coupled to the order parameters with the action: 

^aux. = I dtd''^ (V2s- - ^^s- s + + - -^q+qA (85) 

J V 87r^9s STT^gq J 

Variations 5S'aux. = imply 

g± = 87r\0±, s = M^g,4>lV2 (86) 

Plugging this back into the action, one finds that the interaction is recovered if 

g, - 3gs/2 = -g (87) 

The effective action for the auxiliary s, q^ fields follows from performing the 
fermionic gaussian integrals over the x fields. Let us pass to euclidean space with 
the usual prescription t — > —it, iS —S. We will refer to the D = d + 1 euclidean 
coordinates as simply x. Then the effective action S^q is defined as 



where S'free is the free action for the x fields. The result is 

1 ^ ^ 1 ^ _\ 1 



Sos = d'^xi -^s- s + ^^q^q- ] - ^TrlogA (89) 

where the operator A is a differential operator that depends on s,q^. For constant 
s and g^, the TrlogA can be computed in D-dimensional euclidean momentum 
space since the derivatives are diagonal: {p\d\p) = ip{p\p). For constant fields it 
is meaningful to define the effective potential Vcs = Scs/V^^'^ where V^^'^ is the 
/^-dimensional volume. At finite temperature for instance V^^^ = V (3 where V 
is the usual d dimensional volume and (3 the inverse temperature. Using {p\p) = 
V^^y{27r)^, one obtains 

1 ^ ^ 1 + _ 1 f d^p 



^eff = ^+ -o hdw Triog A(p) (90) 
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In the basis (X| , Xi ) xX) anti-symmetric matrix A is the following: 



A{p) 



/ p2 

— — + Sz 



+ m — Sz 


-x/2s+ 







(91) 



Using Tr log A = log Det A, one finds 
1 ^ 



eflf 



-s ■ s 



s ■ s 



) (92) 



The gap equations follow from setting the variation of V^s with respect to s and 
separately equal to zero. The result is the following for m = 0: 



± 



(93) 



(27r)^p4^g+g- - 

There is no simultaneous solution with both and s non-zero unless one fine tunes 
to the 50(5) invariant point Qs = —Qq- In fact, there is no true competition be- 
tween AF and SC order in these equations, and the distinction between Qs and Qq is 
somewhat fictitious, since for one sign of the coupling there are AF solutions and no 
SC solutions, and visa-versa if the sign is flipped. We thus consider solutions with 
either pure SC order (s = 0) or pure magnetic order (g = 0). To further clarify the 
structure of the gap equation in the sequel, let us separate the spatial and temporal 
parts of the D-dimensional momentum vector p a,s p = (ci;,k). (This notation is 
different from that of section II where there k was a physical momentum; here and 
henceforth it is relative to the Fermi surface.) We also restore the mass m. For pure 
SC order one then obtains the gap equation: 

duj rf-^k 1 



(27r)^+i (cu2 



2^2 



(94) 
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where = q = q. For pure magnetic order one instead has 

_ 2 r dud'^k 1 

1 - -871 Qs J (^2 + k2 + ^2)2 _ ^2 (9^) 

It is important to note the asymmetry in the signs of the above SC versus magnetic 
gap equations, which is ultimately traced to the signs in the 50(5) invariant fl68l) . 
Thus, although the model has an 5*0(5) symmetry that rotates the SC and magnetic 
order parameters, the gap equations are not invariant under the exchange of s and 
q. This implies that the AF gap is related by symmetry to a conventional s-wave 
gap obtained when one flips the sign of the coupling, described in the next section. 

Let us consider now the implications of Goldstone's theorem. A non-zero vac- 
uum expectation value in any one direction of the 5- vector $ preserves an 50(4) 
subgroup of 50(5). Since dim(50(5)) — dim(50(4)) = 4, there are potentially 4 
Goldstone bosons. One of these is associated to the U{1), and when coupled to the 
electromagnetic field is eaten up by the Anderson-Higgs mechanism. This leaves a 
spin-1 triplet of Goldstone bosons associated with the 3-vector 0, and these are the 
closest thing to spinous in our model. The effective theory for these modes follows 
from the above effective potential fl92l) for s. For non-constant fields it contains the 
kinetic energy term d^s- d^s. 

IX. CONVENTIONAL S-WAVE SUPERCONDUCTIVITY IN 3d FOR AT- 
TRACTIVE INTERACTIONS 

It is well understood that a gap that spontaneously breaks the f/(l) symmetry in 
the presence of electromagnetic gauge potential automatically has the characteristic 
electromagnetic properties of a superconductor, i.e. Meissner effect, etc. This pre- 
cursor to the Higgs mechanism can be understood in the original Ginzburg-Landau 



theory 3l|]. A clear explanation of this can be found in Weinberg's book 
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We are primarily interested in positive coupling g since this corresponds to re- 
pulsive interactions. (In the next section we describe how this should correspond to 
the positive U Hubbard model.) However as a check of our formalism thus far, let 
us consider attractive interactions, i.e. negative g, or negative U Hubbard, which 
should reveal the usual s-wave SC instability of the BCS theory. As our analysis 
shows, symplectic fermions give a proper field theoretic description of conventional 
s-wave superconductivity which has not been considered before. 

As expected there is no s-wave superconductivity for positive g since the coupling 
gq is negative and the gap equation of the last section has no solutions. For pure SC 
at negative g, equation ( IHTl) then implies gq = —g is positive. Specializing to d = 3 
one should recover some of the basic features of the BCS theory. The gap equation 
reads 

- = log(l + At/q') (96) 

9q 

As the cut-off goes to infinity, the divergence can be absorbed into the coupling by 
defining l/gq{Ac) = log Ac. This gives the 1-loop correction to the RG beta function: 
l3{gq) = —dgq/ log Ac = g^ and is consistent with the results in Q. A more 
complete RG prescription will be described in section XIII where 1-loop diagrams 
are computed. Note also that positive gq is marginally relevant, i.e. gq increases as 
the energy is lowered. Since q has dimension two, let us define a gap A = with 
units of energy. When the cut-off is large compared to A the solution to the gap 
equation is approximately 

A = ^ = Aee-i/^' (97) 

which is characteristic of the BCS theory, i.e. A vanishes as gq goes to zero, and 
saturates at the cut-off when gq goes to infinity. The zero temperature gap as a 
function of gq behaves as in Figure [61 

Let us now restore the infra-red regulator mass m and confirm the idea that it can 
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FIG. 6: s-wave gap as a function of the coupling. 

be viewed as proportional to the temperature, as long as the latter is close enough 
to zero. As described in section VII, m = aT where a = 1x^^2/3 in 3d. Restoring 
the mass, the gap equation now reads in 3d: 

- = - log ^2 + ^4 J + ^ [tan-i(mVg) - tan-i((m^ + kD/q)] (98) 

When m is too large, the RHS can go through zero and change sign, at which point 
the solution is lost. Numerical solutions of the gap A as a function of m is shown 
for several values of g in Figure [71 

Let mc{g) denote the value of m where the gap disappears. Then, one expects 
the critical temperature to be given by a formula of the form 

T. = ^ = ^Ao (99) 

a a 

where Aq is the gap at m = and Cg is a constant of order unity. This is similar to 
the relation in the BCS theory where Tc ~ Ao/1.76. The constant Cg = mc/Ao can 
be estimated numerically and shows a weak dependence on ^f, as Table I shows. 
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FIG. 7: s-wave gap as a function of the mass m. 
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X. ANTI-FERROMAGNETIC PHASE (MOTT-HUBBARD INSULATOR) 



A. Symplectic fermions from the Hubbard model at half-filling 

As stated in the introduction, one cannot derive our model by taking a direct 
scaling limit of a lattice fermion model like the Hubbard model; this should be clear 
from the fact that they have different symmetries, i.e. 5*0(4) versus 5*0(5). The 
main reason for this is that one needs to take special care of the Fermi surface in 
taking such a continuum limit. Nevertheless, we can further motivate our model by 
comparison with lattice models as follows. Since our model has no explicit lattice, the 
meaning of the s-s^O phase can only be understood by c omp aring with low energy. 



continuum descriptions of magnetism. It is well-known |32l. l35l| that excitations above 
the anti-ferromagnetic Neel state of the Heisenberg model are described by the non- 
linear 0(3) sigma- model. The basic facts about lattice models we need for this 
discussion are collected in Appendix B. Since the n field corresponds to our 0(3) 
vector of order parameters (j) in eq. (l65l) . let us use that notation. 

At half-filling and strong coupling, the Hubbard model can be mapped to the 
Heisenberg model and the order parameter is constrained to have fixed length 
■ = constant. In our model, the magnetic order parameter can be expressed as 

= -^x'^x^ (100) 

where a are the Pauli matrices. A fixed length constraint on (p is equivalent to a 
similar constraint on the x's. Using 

^ij ■ ^ki = 26ii6jk - SijSki (101) 

one can show 

0-0 = -^(x^X+)' (102) 
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where x = Tl,a=^ iXaXa- Therefore if one imposes 

Xf + X'iXt = ^hAc (103) 

for some constant h, then ■ = 3/i^A^/2. (The i is consistent with the propagators 
in Minkowski space, eq. (??)•) 

At half-fiUing there is a simple argument that leads to symplectic fermions. Im- 
posing the above half-filling constraints in the lagrangian: 

d(j) ■ d(f) = AihKc {dx~ dx^^ + irrelevant operators (104) 

where the additional irrelevant operators can be expressed in terms of dimension 4 
current-current interactions J^J'^. Thus, due to the second-order nature of the non- 
linear sigma model, one obtains a fermionic theory that is second order in space and 
time derivatives. Furthermore, let us recall now that we showed in sections II and 
III how symplectic fermions also arise when one expands around a circular Fermi 
surface. Since for lattice fermions the nearly circular Fermi surface occurs inside the 
half-filling diamond, this shows how symplectic fermions can actually extend below 
half filling. Above half-filling (electron rather than hole doping) the same formalism 
applies where now one expands around the circular Fermi surface centered on the 
node k = (vr, vr). 

Relaxing the constraint (11031) moves us away from half filling. A meaningful 
measure of the degree of filling that we will utilize in the sequel is the following. 
Specializing the formula (IHT]) to 2d in euclidean space and with cut-offs, one finds to 
zero-th order: 

1 , , , f^^ d^p 1 1 f A\ , , 

,105, 

Thus, increasing A/Ac corresponds to lowering the density below half- filling, so that 
/i is a measure of hole doping. We will return to this point when we analyze the 
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phase diagram as a function of doping in section XIV, and will also include 1-loop 
corrections. 

The coupling g should be proportional to the Hubbard coupling U (See Appendix 
B). It can be estimated below half-filling near the circular Fermi surface. From the 
lattice kinetic energy eq. (I14p . e(k) ta^k^ for small which gives a Fermi velocity 
vp = 2ta?kF- Since g has inverse length units, U = vpg has units of energy. This 
gives 

^ 2t{kpa) a ^ ^ 

where a is the lattice spacing. 



B. Solutions of the gap equation 

We now study the solutions of the pure AF order gap equation (!95l) at m = 0, 
which should correspond to zero temperature. Because of the two additional minus 
signs in the AF gap equation in comparison with the attractive SC one studied in the 
last section, the interpretation of its solutions is somewhat subtle, and perhaps the 
most delicate point in this whole paper. First of all, the gap equation definitely has 
solutions for gs positive due to the two compensating minus signs. However because 
of the pole at = we require that s > A^ where Ac is the upper momentum cut- 
off. A proper understanding of the solutions requires the introduction of both UV 
and IR cut-offs and to implement the RG directly in 2d. The precise description of 
our RG prescription is postponed until we calculate some explicit Feynman diagrams 
and is described in section XIII. 

When gq = the equation fl87|) gives the relation between g and gs and they are 
both positive. However the resulting relation is purely classical and specific to the 
auxiliary field construction. Instead we fix the relation between gg and g by requiring 
consistency with the perturbative RG. To lowest order this requires gg = 2g, in order 
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for eq. (11101) to be consistent with eq. fll4ip . The gap equation now reads 

1 r^^ »2 



Let us define 



The result of doing the integral is 



S^Al (108) 



^ ^^logf|i±j)-tan-M/<5j (109) 



Since the log-term is positive in the above equation, there are solutions for 6s > 1- 
In fact, when the coupling is lowered, the log term dominates, and the solution 
approaches 6s = and remains there for arbitrarily small g. This behavior is not 
physically sensible since the gap should vanish as the coupling g vanishes. To resolve 
this puzzle, first note that when s = the gap equation has IR divergences that need 
to be regulated. Let us therefore introduce a low-energy cut-off A. Setting s = in 
eq. (I107p . the result can be written as 

^ A A, 

The left hand side is precisely an expansion of the running coupling g{A.), since by 
eq. (]137p . 1/5'(A) ^ (1 + 8g/ A)/g. Denoting the solution as gAF, 

gAF{A) = Y (111) 

Our interpretation of this value of the coupling g^r is that s should vanish for g < g^r 
since s = is a consistent solution at g = gAF- 

The value of the coupling gAF is closely related, but not the same as the low- 
energy fixed point value of g. (See in section XIII.) Namely, at 1-loop the fixed 
point is at g^ = g^K where g^ = 1/8, thus the low energy fixed point occurs at 
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9* ~ T'dAF- As we will see, the low energy quantum critical point at (7* is a second- 
order continuous phase transition that terminates the super-conducting phase on the 
over-doped side. This leads us to propose that the termination point g^p of the AF 
phase is a first-order transition, i.e. the gap drops discontinuously to zero. 

Let us return now to solutions of the gap equation with s 7^ 0. The solution in 
the two asymptotic limits g/ h^c 0, 00 have simple expressions. It is useful to use 
the identity | log(5s + l)/(^s "~ 1) = tanh~^ l/^s- When (7 — 00, 5^ is large. Using 
tanh"^ l/5s — tan~^ 1/5 2/(35^), one obtains 

{glK-^) (112) 

As g decreases, Sg saturates to 1 since the argument of the tanh~^ must be less 
than 1 otherwise the solution is complex. Where the solution starts to flatten out can 
be approximated by extrapolating eq. flll2p down to 5^ = 1, i.e. around g/ = 3/8. 
Using tanh^^ l/6s ~ | log(2/ (6s — 1)) when 6s ~ 1, one obtains 

6s ^1 + le^^^l'"', [g/Kc < 3/8) (113) 

This behavior is shown in Figure [H] where we have incorporated the drop to zero at 
9af- 



XI. D-WAVE GAP EQUATION IN 2d 

In the last section we derived gap equations in the approximation that the gap 
has no momentum dependence. This is essentially equivalent to assuming the lowest 
order scattering of pairs is momentum independent. In Appendix A we show how to 
incorporate momentum dependent scattering and derive the following form of gap 
equation: 
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FIG. 8: Solutions to the anti-ferromagnetic gap equation. 



where the kernel G is given by a Green function related to the scattering of pairs 
with momenta ±k and ±k'. We have set the mass m = 0; it can be restored by 
uj'^ — s> cu^+m^. (As explained in Appendix A, q{k) is not simply the Fourier transform 
of q{x).) 

For our model the kernel will be computed to 1-loop in the next section. In this 
section we analyze the orbital properties of the gap in 2d in a model independent 
way based only on the structure of the gap equation and its symmetries. Similar 
arguments apply to a BCS type of gap equation. 

Let us assume that G is symmetric, G{k,k') = G(k',\<.). If G is also rotationally 
invariant, then its angular dependence arises only through the dependence on k-k' = 
kk' cos{6 — 6'). The kernel and gap can thus be expanded as follows: 



oo 



G(k,k') = ^G'f(A;,A;')cos£(e-^') 



oo 




(115) 
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Substituting the above expansions into the gap equation one sees that due to the 
non-hnearity different I can mix. For simphcity consider a single channel, i.e. assume 
G has only one term in its expansion at fixed I. The angular integral / dO' can be 
performed and turns out to be independent of t 

The result is 



*(^) = -77A^ /"^^ r dk'k'G,{k,k')^—[l- 
^^y J^oo Jo qe{k') y 



(117) 

It's important to note that although G(k, k') varies in sign due to the oscillating 
cosine, the sign of Ge{k,k') is meaningful and determines whether the i channel is 
attractive or repulsive. Negative Ge corresponds to an attractive channel. Further- 
more, any £ = 0, 1, 2, .. is in principle allowed. 

The channel i = 2 can arise rather naturally from a term in the kernel of the 
form — 2(72(k ■ k')^ = —g2k'^k'^{l + cos 26') for g2 a constant, which gives rise to both 
£ = 0,2 with the same sign. As we will show in the next section, for our model £ = 2 
is the first attractive channel. In particular G2 has the form 

G2ik,k') = -87r^g2k^k'^ (118) 

with g2 a positive constant which we will calculate in the next section. This leads to 
a solution of the pure d-wave gap equation of the form 

q{k)=6lk'cos2e = 6l{kl-kl) (119) 

where 6q is a constant satisfying the integral equation: 

54 = 2(72 / dujde i 1 - u^ + k^ ^ I ^^20) 

/(^2 + ^2)2 + 54^4 
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The dependence on k for £ = 2 is of the same form as a particular hnear combination 
of £ = 2 spherical harmonics in 3d, thus we refer to it as dx2-y2, or simply d-wave, 
as in the literature; i = 0,1 can be referred to as s and p wave. 

The above gap equation has some interesting properties, in particular, 6q = 
when g2 is too small. Since this kind of gap equation must be regularized in the UV, 
we are led to define 

^2 = 1 (121) 
where ^2 is dimensionless. To estimate the lowest value of ^2 with non-zero gap, the 
integrand in the above equation can be expanded in powers of Sg. Keeping terms of 



order 6^, one finds that for 6q small it behaves as 



a' V «fl'2 



1/4 

(122) 



where 

a = (8 + 7r -81og2)/12, a' = 3(104 + Stt - 128 log2)/384 (123) 

Thus 

5g = for ^2 < - ~ 2.15 (124) 

a 

For ^2 large, one finds 

^ {292)'" (125) 

However as we will explain in section XIV, there will also be an upper threshold that 
comes about when the RG is properly implemented, and this leads to a SC dome. 

XII. 1-LOOP SCATTERING AND THE GAP EQUATION KERNEL 
A. Feynman diagrams 



In this section we derive the kernel for our model and show that it has an attractive 
d-wave channel. Let G^^\pi,p2^'P'i^'Pi) denote the 4-particle vertex function with 
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the overall (211)^ 5^^\p) removed. (Apart from this overall factor, G^^^ differs from 
r*^^'' in Appendix A by an overall sign since the Feynman rule vertex is defined to 
be negative in euclidean space.) The arrows in the figure indicate how the flavor 
indices a,(3 = 1,..,N are contracted. They also indicate the flow of charge since 
the interaction is proportional to (XlaXoXa)^- This structure is shown in Figure [9], 
which indicates the interaction of the euclidean space Feynman rules. 




FIG. 9: Interaction vertex 

The 1-loop contributions to G^^^ are displayed in Figure [10] and consist of 3 sep- 
arate channels which differ in how the a, (3 flavor indices are contracted and also in 
their momentum dependence. It is interesting to carry out this part of the calcula- 
tion for arbitrary A^. The Sp{2N) group theory factors are 2 — A^, 1, 1 respectively 
for the three diagrams, where the — A^ dependence comes about from the closed loop 
and a fermionic minus sign. 




FIG. 10: 1-loop Feynman diagrams with group theory factors 
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To this order one thus has 



G^'\puP2,P3,P,) = ^n'g - 2{4n'gr [(2 - N)f{pl,) + f{pl,) + f{p 

where Pij = Pi+ Pj and / is the function: 

d^i 1 



(126) 



fip' 



(27r)« [P + m^][{£ + py + 



(127) 



d 



(27r)« 



dx- 



1 



+ x(l — x)p'^ + m^]" 



Since we are interested in low energies, it is meaningful to expand the integrand 
in powers of p^ ji^ . At zero temperature, we can also set m = 0. One finds 



f{p^ 



d^i \ ( ^ Ip'^ _ I p^ 



:i28i 



(27r)^ 3£2 ' 10£4 ' • 

As explained in Appendix A, the gap equation kernel is obtained by specializing to 
pairs of opposite momenta. We thus fix 



Pi = -P2 =P, P3 = -Pa = P 



(129) 



There are also 2 inequivalent ways in which to contract the vertex with external legs 
which leads us to define 



G{p,p') = G^'\p, -p,p', -p') + {p' ^ -p') 



(130) 



The momentum independent part is the following: 



^(^,^01^^=0 = 87r^(7 - 647rV(4 - N) 



d 



(131) 



The second term in the above equation can be absorbed into a redefinition of the 
coupling g. It in fact is just the 1-loop contribution to the RG beta-function and we 
will return to it in section XIII where our RG prescriptions will be fixed. 
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The order term is 



647r4(3 - A^) 



d 



D, 



(27r)^£6 



(132) 



According to the derivation in Appendix A, for static gaps, the kernel G(k, k') in 
the gap equation is simply G{p,p') with the time components of p disregarded, i.e. 
p = (0, k). This implies that though ^(k, k') is closely related to the S-matrix, it is 
not identical since the external momenta are not on-shell. The above term is then 
a repulsive correction to the s-wave contribution, which is already repulsive at tree 
level. 

Finally we come to the p^ term: 

327r^(3 - N) 



G(k,k')|,4 



J2\2 



4(k ■ k')') 



d 



(27r)^£8 



Using 



(k ■ k')' = k^k'^ cos\9 - 6') = ek'\cos2{9 - 9') + l)/2 



(133) 



(134) 



one sees that this term gives both s and d-wave contributions. Whether the channel 
is attractive depends on the sign of the factor 3 — A^. For N = 2, the above 1-loop 
s-wave correction is attractive but since the leading tree-level term is repulsive, it is 
unlikely that it could lead to s-wave pairing. The leading contribution to the d-wave 
term is attractive for < 3. This is rather interesting since this d-wave instability 
is invisible to large methods. In the sequel we will make the approximation of the 
last section and not consider possible mixing of the s and d-wave gaps. 
The d-wave term G2{k, k') is of the form in eq. f lllSp with 
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(3-A^)87rV f 1 
(27r)3 £8 



(135) 



The gap thus has the characteristic d-wave form in eq. (11191) where the constant Sq 
is a solution to eq. (11201) . 
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XIII. RENORMALIZATION GROUP SPECIFICALLY IN 2d 



In the present context it is most appropriate to adopt the Wilsonian effective 
action view of the RG. In this approach there are two energy scales to consider, a 
fixed cut-off scale A^, and a running scale A < Ac- Integrating out high energy modes 
leads to couplings that depend on the running scale A. Wherever possible, one lets 
the cut-off Ac go to infinity. 

For illustration, consider first the contribution of loop modes with A < £ < A^, to 
the momentum-independent part of G in "id: 

G{p, p') = Stt' (^g-{A- N)g' £^ ^) (136) 

The effect of these high energy modes is to modify the coupling according to g{A) = 
g + g'^{A - iV) log A/Ac. The beta function is then -dg/d\ogA = -(4 - A^)^'^, 
in agreement with results injl?!. 18]. An equivalent manner to describe this RG 



prescription is to let J di ^ J^^" di and then set the cut-off Ac — > oo and keep the 
dependence on the renormalization energy scale A. 

Let us now specialize this to 2d. The momentum independent part of G is now 
the following 

G = Svr^ {g - 4(4 - N)gyA) ^ Svr^ ^(A) (137) 

where we have sent the cut-off Ac to oo. This gives —dg/dlogA = —4(4 — N)g'^/A, 
which is consistent with the fact that g has units of energy. The latter is an important 
feature of our model, and in order to deal with it, we define a dimensionless coupling 
9- 

g{A)=AgiA) (138) 

The beta function for ^ is now 

-A^ = ?-4(4-iV)r (139) 
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where the hnear term just reflects the classical dimension 1 of g. There is a fixed 
point at 



1 



9* 



(140) 



4(4 -iV) 

where the above is only approximate due to corrections beyond 1-loop. 
The solution to the RG fiow equation is then 

- A.4(4-'4k-A)go 

The fixed point value is reached irregardless of whether the initial coupling at 
the cut-off Ac is large or small, as long as it is positive. This behavior is sketched in 
Figure [Til 




A 



FIG. 11: RG flow to | at low energies from either strong or weak coupling at short distances. 
The coupling g2 that enters the d-wave gap equation can now be expressed as 



92 



25 A3 



(142) 



where we have safely taken Ac to oo in eq. 01351) and set N = 2. 
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XIV. SCALING AND GLOBAL FEATURES OF THE PHASE DIAGRAM 



A. RG scaling 

For the purposes of comparison with real compounds, our model has only two 
independent parameters. First is the strength of the coupling at the short distance 
scale Ac- Second is the overall scale A^. The complete phase diagram falls into 
place when one takes into account the RG scaling. Throughout this section it is 
implicit that we are only working to 1-loop and all equalities should be taken as 
approximations. 

The structure of the phase diagram is most transparent and the universal features 
more clearly revealed if one works with the dimensionless inverse couplings: 

111 . , 

X = xo = —, x^ = — (143) 

9 9o 9* 

Throughout this section we fix = 2 where x* = 8 is the fixed point value. The 
scale A of g = Kg can now be expressed in a simple way in terms of x using eq. 
(HH]): 

A ^ 

Ac X* - Xo 

The variable x should not be confused with the conventional doping variable in the 
literature, although it is closely related; the hole doping is the variable h below. The 
above linear relation is specific to 2d and is crucial to understanding the overall struc- 
ture of the phase diagram. Since the line A/A^ represents the scale of the coupling, 
below this line energies are comparable to the coupling and the non-Fermi liquid 
properties begin to reveal themselves. We have labeled this region the pseudogap 
as in the literature and will return to it in section XVI. Note that when x = xq, 
A/Ac = 1. 

As we defined them, both the AF and SC gaps 6s and 6q are in units of the cut-off 
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Ac. Since we have performed an RG transformation of the couphng from Ac to A, 
we must also scale the gaps as follows: 

It will also be convenient to define the parameter 7 which measures the distance to 
the fixed point at short distances: 

7 = ^1^^ (146) 

where < 7 < 1. The AF gap equation fllOQp should now be solved for S'^i^x) where 
g is replaced by 

f = - f (147) 

Ac X \X^, — 2^0 / 

There are two cases to consider depending on whether the coupling ^ is strong or 
weak at short distances, i.e. whether ^0 is above or below the fixed point value 
corresponding to the two curves in Figure [11], which we will refer to as Type A and 
B. 

Type A: strong coupling at short distances 

Here we assume 'go > 'g* = 1/8, which implies xo = {1 — 7)x=i.. The scale as a 
function of x now has negative slope 

A . _i (± _ 1) (148) 

Ac 7 V^;* / 

and since the above ratio is by definition less than 1, we have Xq < x < x^,. From 
eq. f illip the first-order transition that terminates the AF phase corresponds to 
dAF = ?*Ac/A which translates to 

XAF = (149) 

1 + 7 
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The value xaf is always to the right of xq since 7^ > 0. 

Let us turn now to the SC phase. The d-wave gap equation (11201) depends on (72, 
which by eqns. (11211 11421) can be expressed as 

= 25x2(1 -x/x.)3 

The d-wave gap equation only has solutions if g2 is positive, which requires x < x^.. 
Furthermore, ^2 niust be above the threshold ^2 > 1/a by eq. (I124p . Thus the d-wave 
SC gap is non-zero in the range 

5g 7^ for xi < X < X* (151) 

where the lower threshold is a solution to 

4a7^ = 25x^(1 - xi/x^f (152) 

(The constant a ^ .466 is defined in eq. (I123p .) It should be emphasized that, unlike 
the values of Xq, xaf and x^, which are determined by the low energy fixed point, the 
value Xi is not universal, i.e. not predicted by the RG itself. Since the solution to 
the above equation is close to x*, one can estimate 

4a ^ ^''^ 



xi^x, I 1-7(^^J l^(l-0.l7)x, (153) 

In principle it is mathematically possible to have a non-zero solution to the d-wave 
gap for X sufficiently small, but this will tend to be small and inside the AF phase. 

The transition point xaf is outside the SC dome in the approximation we have 
made, since xaf > Xi would require 7 > 1. These features, and the geometrical 
relationships between the various transition points, are shown in Figure [121 which is 
not drawn to scale since the latter depends on 7. 

In Figure [13] we display numerical solutions to the AF and d-wave SC gap equa- 
tions for 7 = 1. Rough comparison with experimental results suggests the SC dome 
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FIG. 12: Global phase diagram for Type A. The slope of the dashed line is ——zr- The 
various transition points are related geometrically by xq = (1 — 7)x=i,, xaf = + 7), 

and xi (1 — .17)x=k, where x* = 8. The dashed line represents eq. (jl44p . 

as we have calculated it is too narrow, which would imply that we overestimated 
Xi- However as stated above, the value of Xi is less universal than xq^xaf and x^,, 
and could easily change by improving our approximations, for example taking into 
account s and d wave mixing, or incorporating other effects we have neglected, such 
as interplane coupling and disorder. The important point is that there is both an 
onset and termination of the SC phase, and the termination point on the overdoped 
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side is what is universal since it corresponds precisely to a quantum critical point. 
Furthermore, since xaf < 2^1 , AF and d-wave SC do not appear to compete so that 
the SC is robust. 




X 



3.0 4.0 7.1 8.0 



FIG. 13: Numerical solutions of the AF and d-wave SC gap equations for 7 = 1 (Type A). 

Type B: weak coupling at short distances 

In this case ^0 < 9* at short distances, and 7 = (xq — where in principle 

< 7 < 00. Now the slope of A versus x is positive: 

K I f X \ , , 

- = - - - 1 (154) 



Ac 7 \x^ 
and x^ < X < xq. The AF transition point is at 



Xaf = (155) 
1-7 
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Thus < 7 < 1 otherwise xaf goes from oo to — oo. 

There are non-zero solutions to the d-wave gap equation in the range: 

(5g 7^ for < X < X2 (156) 

where X2 is a solution to eq. (11521) with xi — »■ X2 and 7 — >■ —7. An approximation 
of the form (1153^ also holds with 7 —7. With the range of parameters of our 
model, Xaf lies outside the SC dome, so that the SC dome is completely inside the 
AF phase. These features are shown in Figure O 



B. Doping: 1-loop corrections and optimal doping 

The hole-doping variable h defined in section X takes the simple form in terms of 
the X variables: 

h{x) = 1 ( ^) (157) 



TT^ \X^ — Xq 

Again, this linear dependence on x is characteristic to 2d. As argued in section X, 
since increasing A/Ac decreases the density to below half filling, /i is a measure of 
doping, where x = Xq should correspond to half-filling. Thus apart from an overall 
scale of 1 /tt^ and a shift of the origin, our previous phase diagrams in terms of x are 
effectively in terms of the doping h. 

Since both the AF and SC termination points xaf and are closely related to 
the RG fixed point, the ratio of doping at these two values should be universal. Let 
Haf = h{xAF) and /i* = /i(x*). Expressed in terms of 7 one has 



1 7 1 
hAF = , /i* = — J (15^ 

TT^ 1 + 7 TT^ 



which gives the ratio Haf/K = 7/(1 + 7). 

The above formula for h is from eq. (11051) where only the zero-th order contri- 
bution was calculated. We can easily include the first order correction in g, which 
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FIG. 14: Global phase diagram for Type B. The values shown are related as follows: 
xq = (1 + 7)3;*, XAF = — 7), and X2 ~ (1 + .l7)x=i,, where = 8. The dashed line 

has slope l/{'yx^). 

is the 1-loop self-energy diagram shown in Figure [151 This leads to the following 
correction: 

There is an important fermionic minus sign in the above expression coming from the 
fermion loop. Expressing the result of doing the above integral in terms of x using 
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eq. (I144p one obtains 



h{x) 



X — Xq 



1 + 



4 / X — Xq 



71 \X^ — Xq ^ 

The corrections to the formula (11581) are then 

1 7(2 + 7) 



h 



AF 



(160) 



(161) 



7r2 2(1 + 7) ' 27r2 
The optimal doping fraction can be estimated as follows. Recall that for Type A, 
the SC phase occurs for xi < x < x^. Evaluating h at the lower limit Xi one finds a 
weak dependence on 7: 

h{xi) 



1 /1.3- 0.097 
~ \ 1 - O.I7 



One thus concludes that optimal doping occurs in the tight range: 



1.3 _ _ 3 

— 2 — -13 < /^optimal < .15 — 



vr 



27r2 



This leads to 



^ 7(2 + 7) 

/^optimal 3(1 + 7) 



(162) 



(163) 



(164) 



For 7 = 1/2 this gives Haf ~ -04. 



+ 



A 12^1^ 



FIG. 15: 1-loop correction to the propagator 



XV. DETAILS ON THE AF AND D-WAVE SC GAPS AND CRITICAL 
TEMPERATURES 

In this section we provide detailed numerical solutions to the gap equations for a 
variety of 7 and also reintroduce the mass to incorporate a temperature. 
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A. AF gap 



Numerical solutions to the AF gap equation (llOOp expressed in terms of x for 
7 = 1/2 are shown in Figure [TBI 




Xo 



XAF 



FIG. 16: Numerical solutions to the AF gap for 7 = 1/2. xq = 4 and xaf = 5.3. 

To study finite temperature, we introduce a mass m = aT where a = tt^/^/ \/6 ~ 
1.7 was introduced in section VII. The gap equation becomes 

^ = 1 I ^52-m2tanh-i , ^ - 6^, + mHa.n~^ , ^ | (165) 

where m = m/Ac- Since we have RG scaled the gaps, in order to make comparisons 
on the same scale we need to also define a scaled mass: 

, A , . 

m = —m (166) 
Ac 

Numerical study of eq. (11651) shows that the solution S'^ vanishes when m' is 
too large, consistent with identification of m with temperature, as in conventional 
superconductivity discussed in section IX. Let m'^ denote the value of m' where 
the solution 6^ vanishes. This leads us to define a Neel temperature T/v where the 
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gap vanishes. As in the BCS theory, Tjv is expected to be proportional to the zero 
temperature gap. Let us define then 



and in these equations 6'^ is the zero temperature gap. Inspection of the finite tem- 
perature gap equation f ll65p shows that the solution should disappear when the argu- 
ment of the square-root is negative, i.e. when rh ^ 6s, and caf ~ 1- This is a delicate 
limit, however we have verified that this is approximately correct numerically. 

B. SC gap 

Numerical solutions to the d-wave gap equation fll20p expressed in terms of x 
are shown in Figure [17] for 7 = i, |, |, 1. One sees that the peak value of the gap 
(5g ~ 0.11 is not very sensitive to 7. 

Let us now introduce a temperature T by letting cj^ u'^ + m'^ in the gap equation 
fll20p where as before m = aT. The behavior of the gap as a function of T is shown 
in Figure [18] for 7 = i, i, 1. Let m'^ denote the value of m where the gap vanishes. 
One sees from the figure that here there is a somewhat stronger dependence of 
on the inverse coupling x than for the AF gap. 

The critical temperature Tc can now be defined as 




(167) 



a 



a 



where 




(168) 




a 



a 



(169) 



where 




(170) 
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7.0 



= 8 



FIG. 17: Numerical solutions of the d-wave gap equation for 7 = |, ^, |, 1. 

is a dimensionless constant. In this final formula, we have restored fundamental 
constants and the Fermi velocity vp. Above, 6^ is the zero temperature gap, which 
depends on 7, x as does m'^. 

Let Xopt denote the value of x with the largest gap. The constant csc for various 
7 at the corresponding Xopt are shown in Table 2. One sees that csc is of order 1 
and has only a weak dependence on the couplings. 
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m! Ik, 

0.1 



FIG. 18: d-wave gap as a function of ra for 7 = 1/4, 1/2, 1. 



7 




<^g(a;opt) 




csc = m'^/S'gAc 


1 


7.5 


0.118 


0.06 


0.51 


3/4 


7.6 


0.116 


0.06 


0.52 


1/2 


7.8 


0.114 


0.07 


0.61 


1/4 


7.9 


0.113 


0.08 


0.71 


1/8 


7.94 


0.113 


0.06 


0.53 


1/16 


7.97 


0.113 


0.06 


0.53 



We can estimate now the critical temperature at optimal doping. The above table 
shows that csc depends relatively weakly on 7. Furthermore, 6'^ at optimal doping 
is also relatively constant ~ 0.11. Thus the scale of the critical temperature fll69p 



70 



is set primarily by vp and Ac, i.e. it is only weakly dependent on the coupling g or 
the Hubbard couplings by eq. (11061) : this is the beauty of having a low energy fixed 
point. The scale Ac = 1/a is an inverse length a, which should be on the order of 
the lattice spacing. 

To give a rough estimate of Tc, let us take the lattice spacing to be that of the 
CuOi square lattice a = 3.8A. For vp we use the universal nodal Fermi velocity 



from[20|, which we estimated to be vp ~ 210km/ s in section III for LSCO. The 
most unknown quantity is a that sets the relation between temperature and mass; 
let use the estimate of a = from section VII. With the average values 

6g = 0.11, csc = 0.6, from the above table, this gives Tc ~ 140-?^, which is quite 
reasonable considering all of the order 1 constants we have approximated. Note that 
this result relies on the relatively small value of 6'^ we found from numerical solutions 
of the d-wave gap equation; a value 6'^ of order 1 would give Tc higher by an order 
of magnitude. 

A useful form of the above equation for Tc is 

Tc = csc — - 650K (171) 
a 

where a is in A and vp in ev — A, and we have set a to its estimated value a = 1.7. 
For the range of csc shown in the above table, 120K <Tc < 160i^. The maximum 
Tc occurs around 7 = 1/4. Since this Tc is on the high side, this is perhaps because 
we underestimated a in our simplified inclusion of temperature. Another possibility 
is that Ac should instead be set by the average separation of holes. At doping /i = .15 
this increases a by a factor of about 2.6 leading to AQK < Tc < Q2K. 

The above formula gives some hints on how to increase T^. shorten the lattice 
spacing, increase the Fermi velocity by somehow modifying the effective electron 
mass m*, or tuning the material to 7 = 1/4, which would require screening the 
Coulomb potential at short distances. In particular Tc should increase with pressure 
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if the main effect of higher pressure is to reduce the lattice spacing. 



XVI. THE PSEUDOGAP REGION 

It is often suggested in the hterature that the pseudogap is the key to understand- 
ing high Tc since the AF and SC order condense out of it, and there has been much 



speculation about its nature. (For a review of the experimental data see 33|].) In this 
section we briefly discuss the insights on the pseudogap furnished by our model thus 
far. 

As stressed throughout this paper, the essential ingredient is the non-Fermi liquid 
at higher temperatures before condensation. The pseudogap is the region where the 
energy scales are such that this non-Fermi liquid behavior is visible, i.e. there are 
strong correlations of the electrons described by our fields but the energy is not 
low enough for them to condense into an AF or SC ordered state. The relevant 
energy scale is A of the coupling g = Ag, thus it is natural to define a pseudogap 
temperature Tpg-. 

Tp, = A = xg{A) = -(i-^)a, (172) 

where we have used eq. (11471) . Thus, as defined, Tpg is proportional the coupling g 
and below this scale the order g corrections are relatively large. Tpg corresponds to 
the dashed line in the figures of section XIV, and this line is seen in measurements 
of various properties on the underdoped side. Our model predicts that the Tpg line 
runs between the top of the AF gap to precisely the quantum critical point which 
terminates the SC phase on the overdoped side, and the experiments support this 33| . 

Since the dashed line in the figures simply represents the relevant energy scale 
where the interactions become strong, it does not represent a phase transition. There 
are thus no spontaneously broken SU{2) or U{1) symmetries just below Tpg. The 
field with non-zero expectation value is the SU(2) ® f/(l) invariant {x~X'^) which by 
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eq. fllOSp is related to A/Ac- The spectrum still consists of particles and holes with 
charge ibe, just strongly coupled. Any experiment performed at energies comparable 
to the scale of the coupling g is sensitive to the quantum corrections in our theory, 
thus there wide variety of probes of Tpg. 

A possible concrete manifestation of the above pseudogap scale is a dynamically 
generated mass m at zero temperature, since such a mass is necessarily proportional 
to A. At finite temperature this mass will typically lead to thermal suppression 
factors e~'^l'^ at low temperatures. This idea will be investigated in forthcoming 
work. 

XVII. CONCLUDING REMARKS 

To summarize, we have constructed a 2d relativistic model with 4-fermion in- 
teractions and analyzed in detail many of its properties which revealed a striking 
resemblance with the known features of high Tc superconductivity in the cuprates. 
Apart from the overall scale Ac, the phase diagram can be calculated as a function 
of a single parameter 7, and gives a very good first draft of what is observed in 
experiments. Since these features of the model were outlined in the introduction, we 
conclude by summarizing the main theoretical aspects that are new and responsible 
for the model's properties. 

• In the approach we have taken for expanding around a circular Fermi surface, there 
is essentially a unique non-Fermi liquid in Id for spin ^ electrons. The requirements 
of a local, rotationally invariant quantum field theory in a sense make this theory 
inevitable. The rotational invariance is not in conflict with the existence of a lattice 
if the wavelengths at low energies are long compared with the lattice spacing, in fact 
it is already known that the non-linear sigma model description of the Heisenberg 
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magnet is rotationally invariant. 



• Our model has a low energy fixed point for the same reason that the 0{M) vector 
models have the Wilson-Fisher fixed point, since the perturbative expansion differs 
from the latter only by some fermionic minus signs. The model is perturbatively 
tractable since the low energy value of the coupling is relatively small ~ 1/8, even 
for arbitrarily strong interactions at short distances. 



ar 



'ermion with in- 



• The kind of theory considered here, namely a Lorentz seal 
teractions, has only recently been considered in this context [l7, [l8| because it was 
previously thought to be inconsistent with the spin-statistics theorem and unitarity. 
In fact, the free version arises as Faddeev- Popov ghosts in gauge theory. As explained 
in 1^, since the theory turns out to be pseudo-hermitian, H'' = CHC, it still defines 
a unitary time evolution and has real eigenvalues. The non-interacting theory is ac- 
tually a perfectly hermitian description of particles and holes near the Fermi surface. 
In this paper we further elaborated on this issue, noting that the operator C simply 
distinguishes between particles and holes. Furthermore, the kinematic constraints 
coming from the expansion around the Fermi surface require eigenstates of H that 
are also eigenstates of C, and for such states H is actually hermitian. Finally, the 
model is consistent with spin-statistics since spin is a flavor symmetry and we are 
still quantizing spin ^ particles as fermions. 

• A small non-zero temperature can be introduced as a relativistic mass term in the 
lagrangian. Thus, in euclidean space, the theory is 3 dimensional and spontaneous 
symmetry breaking is possible, in accordance with the Mermin- Wagner theorem. 

• Our model gives a new quantum field theoretic treatment of conventional s-wave 
superconductivity for attractive interactions in 3d and can be extended down to 2d. 
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• The same model can be motivated from the Hubbard model at half-filling. Thus, 
the model can interpolate between the nearly circular Fermi surface just below and up 
to half-filling. This analysis shows how to vary the hole doping by varying the cut-offs, 
or equivalently the inverse coupling, and leads to a calculable phase diagram. The 
phase diagram has some universal geometrical features that depend on the strength of 
the coupling at short distances. In order to properly understand the phase diagram, 
it was essential to recognize that hole doping was proportional to the inverse coupling 
and to implement certain RG scaling relations. 

• It was necessary to derive a new kind of momentum-dependent gap equation that 
incorporates the scattering of Cooper pairs near the Fermi surface. Only then can 
one see that there is an attractive d-wave channel even if the original interactions 
that led to AF order were repulsive. This d-wave channel is only attractive for < 3, 
so it cannot be understood using large methods. 

• Although the 5*0 (5) symmetry helps to explain the existence of both AF and SC 
order since order parameters for both are present, since the AF gap is s-wave and 
the SC gap is d-wave, they are not simply related by symmetry. In particular the 
d-wave gap equation is second order in the coupling g whereas the AF one is first 
order. What is related by 5*0(5) symmetry to the AF phase is a conventional s-wave 
SC obtained when one fiips the sign of the interaction to make it attractive. 

• The pseudogap is the region where the energy scales are comparable to the scale 
of the dimensionfuU coupling but not low enough for the particles to condense. 

Our theory certainly reproduces all of the qualitative features of high materials 



and we believe also explains the universal nodal Fermi velocity observed in [20|]. On 
the quantitative side, the 1-loop calculation of the optimal doping fraction is in good 
agreement with experimental results, as is our estimate of Tc. It is beyond the scope 
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of this article to present additional detailed comparisons with the many impressive 
experimental results obtained over the past 2 decades. 
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XIX. APPENDIX: DERIVATION OF THE MOMENTUM DEPENDENT 
GAP EQUATION 

In order to consistently incorporate higher order scattering, one can start with the 
so-called IPI effective action S^s with the usual definition as the sum of 1-particle 
irreducible verticesjsl. Throughout this section we work in euclidean space. Since 
we are only interested in gaps that depend on the spatial component k of p = {u, k), 
one should start with a 4-particle term of the form: 

•S'effU-particie = dt ci'^xi ■ ■ ■ (i'^X4 T^^\xi, X2, X3, X4) xf (a^OXj^ (^2)%^ (a^s)^!' (^4) 



where all Xj are at the same time: Xi = (t, Xj). However in order to streamline 
the derivation, and also to more clearly express the result in terms of the usual 
momentum-space Green functions, we treat space and time on equal footing and 
make the above reduction at the end. The 4-particle term in our N = 2 model then 
has the form 



5'eff 1 4-particle = J d Xi ■ ■ ■ d Xi T^'^' {Xi, X2, X3, Xi) {Xl)Xl {x2)X~i i^3)xt i^i) 
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(173) 




where (dp) = d^p/ (2,7:)^ and x{p) is the Fourier transform of xi^)- The function 
r^^^ is a 4-point correlation function of x's and by translational invariance has an 
overall (5-function: 

r^'nPi,P2,P3,P4) = (27r)^5(^)(pi+p2+P3+P4) (-871^9 + ....) (174) 



where we have included the tree-level contribution and .... represents loop corrections. 

To derive the gap equation we follow the auxiliary field method of section VIII 
and introduce momentum dependence in the manner described by Weinberg^, [sj. 
Introduce pair fields q^{piiP2) and the auxiliary action 



'S'aux = j {dpi)---{dpi)T^^\pi,p2,p-i,Pi) -q'^{pi,P2)q {P3,Pi) (175) 

+ <l~{PuP2)X^P3)xtiP4) + q^{PuP2)xj{P3)X^{P4 

The equations of motion for q give q~ = XjXj Q.^ = X^X^ ^'^^ substituting 
back into ^aux recovers the correct quartic interaction of x's- 
We now specialize to Cooper pairs with total momentum zero: 

qiPuP2) = (27r)^5(^)(pi+p2)g(pi) (176) 
Define the kernel G as follows: 

T^'\pi,-Pi,P3,P4) = (2vr)^5(^)(p3 + P4)G(pi,P3) (177) 
Then the auxiliary action becomes 

^aux = J{dp){dp')G{p,p')[-V^''k+{p)q~{p') + q^{p)xt{p')xt(-P') (178) 

+ (i~^{p)xj{p')x'i{-p') 

where V^^^ = (27r)^5*^^)(0) is the D-dimensional volume. 
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The free field kinetic term is the following 

^free = [ (dp) J2 Xaip)ip' + m')xt{-p) (179) 

The gaussian functional integral over the x fields can now be performed to give an 
effective potential Ves = Scs/V^^\ To describe V^s in a compact form it is convenient 
to define 

g±(p) = l{dp')G{p,p')q^{p') (180) 

One then finds 

Kff = - y" {dp){dp') G-\p,p') q+{p')r{p) -\j {dp) Tr log A(p) (181) 

where A{p) is of the same form as in eq. (!9T]) with —Q^ip) and s = 0. Above, 

G^^ is the inverse of G as an integral operator, and not simply 1/G. 

In order to understand how one recovers the constant gap equation, let us pass 
to position space where 

q{xi,X2) = j{dp)e'P<^'-''^q{p) (182) 

For a constant kernel, it can be redefined to G{p,p') = 1, and one sees that 

q{x,x) = J {dp)q{p) = q (183) 

thus a constant q{p) corresponds to a constant q = q{x, x) in position space. Thus, 
where there is no time dependence in the gap, one can let p = [uj, k) and define 

g(k) = Jidk')Gik,k')qik') (184) 

where now G(k, k') is the same kernel but with the time Fourier transform disre- 
garded. In other words, (^(k, k') is simply G{p,p') with p = (0, k) and p' = (0, k'). 
Finally, 5gKfr = gives the gap equation 

r dud'^k' , q(k') , , 

= - / G(k, k') ^^^^^^^ (185) 

where q = q^- After relabeling q q this is the equation (I114p . 
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XX. APPENDIX B: RELEVANT ASPECTS OF LATTICE FERMIONS 



In this appendix we collect some known features of lattice 



"ermion models that 



34 



35|. 



are referred to in the paper. All are contained in the reviews 

Let Fj denote the positions of sites of a 2-dimensional square lattice with lattice 
spacing a. Nearest neighbors and are related by = + a with a G {ai, .., 
where ai = — a3 = (a, 0), a2 = —a^ = (0, a). Introduce fermion operators c^^^a where 
a =t, i represents spin and define the hamiltonian 

H=-tYl (4,aCr+a,.) (186) 
r,a,a 

where 

{Cr,Q,4,^^,} = 5r,r'Sa,a' (187) 

The hamiltonian is hermitian due to = Xl-a- 
Introduce the momentum space expansion 

Cr,a = Xle"'"Ck,a (188) 
k 

Using Xlr^*^ "^ ~ ^k,0) the hamiltonian becomes 

k,a 

where 

Ey^ = —2t {cos k^a + cos kytt) (190) 
Define the local spin operators: 

•^r^ = Cri, S; = -^cl^ Crt, = ^ (^cj^ C^^ - cj^ Cr|) (191) 

The conserved spin SU (2) charges are Q = J2r '^r and satisfy the SU (2) algebra 
with the convention in eq. ( l58i) . 
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The Hubbard interaction is essentially unique up to shifts of the chemical potential 
as a consequence of the Fermi statistics: 

H,^, = UY,{n,^-\){n,^-\) (192) 

r 

where rira = cj^, Cra are local number operators. The interaction can also be written 
as a spin-spin interaction since 

^ ^ 3 3 

Sr- Sr = -T^^r^nrl + -{Ur^ + Uj, ) (193) 

Z 4 

In addition to the above spin SU (2) symmetry, there is another commuting SU (2) 
symm etry so that the largest symmetry of the Hubbard model is SU (2) ® SU (2) = 



S0{4:) 36|. The latter SU{2) is intrinsic to the lattice and does not obviously have a 
continuum limit. 

At half-filling there is one fermion per site, which implies nr^nrj = and rij.^ + 
riri = 1. Thus at half- filling one sees from eq. f ll93p 

Sr-Sr = ^ (194) 

Equating the above with j(j + 1), j = 1/2 one sees that at half- filling the local spin 
operators S form the 2-dimensional spin | representation of SU{2). Away from half 
filling the constraint (^r)^ = constant needs to be relaxed. 



At large U and half-fi^ 
spin I Heisenberg model 



ling, the Hubbard model can be formulated as an effective 
H,s = jJ2 ('S'r. ■ 4 - 1/4) (195) 







+2 



where J = Ar/U. 

In the continuum limit an effective theory for fluctuations above the AF state is 



the non-linear 0(3) sigma model 32]. One starts from a staggered configuration 

Sr = ±nr + ir (196) 
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where ± is for even/odd sublattices. In the continuum hmit the effective theory for 
the n field is second-order in space and time derivatives with lagrangian density 

£ = i {dtn ■ dtn - W ■ Vn) + Aop (197) 

where is constrained due to eq. f ll94p 

= constant (198) 

The topological term has the form C^^p ^ n ■ dn x On. In Id the topological term is 
known to serve an important role in determining the low energy fixed point 32|]. In 
2d the topological term is not known to play an analogous role; in fact it is an RG 
irrelevant operator of dimension 7/2 and we ignore it. 
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